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Abstract 

Within the precise cosmological framework provided by the A-Cold Dark Matter model and 
standard Big Bang nucleosynthesis, the chemical evolution of the pregalactic gas can now be 
followed with accuracy limited only by the uncertainties on the reaction rates. Starting during 
the recombination era, the formation of the first molecules and molecular ions containing hydro- 
gen, deuterium, helium, and lithium was severely hindered by the low density of the expanding 
universe, the intensity of the cosmic radiation field, and the absence of solid catalyzers. Molecu- 
lar hydrogen and deuterated hydrogen, the most abundant species formed in the gas phase prior 
to structure formation, played a fundamental role in the cooling of the gas clouds that gave birth 
to the first stellar generation, contributing to determine the scale of fragmentation. Primordial 
molecules also interacted with the photons of the cosmic background via resonant scattering, 
absorption and emission. In this review we examine the current status of the chemistry of the 
early universe and discuss the most relevant reactions for which uncertainties still exist from 
theory or laboratory experiments. The prospects for detecting spectral distortions or spatial 
anisotropies due to the first atoms and molecules are also addressed. 
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1 PRESENT CONTEXT AND MOTIVATIONS 

We live in a chemically and dynamically evolving universe. The interplay be- 
tween these two domains has shaped the universe in the way we observe it now. 
Dynamics, driven by visible and invisible forces, has governed the transformation 
of minute density perturbations into structures of varying degree of complexity 
(from galaxies to individual stars within them). Chemistry has been instrumen- 
tal for modifying a very simple mixture of nuclei and electrons into complex 
molecules and solid particles, the essential ingredients for cooling and fragmen- 
tation of the interstellar gas. This review is an account of the portion of the 
chemical evolution that spanned a small fraction of the universe's lifetime, be- 
tween the recombination era at redshift z ~ 1000 when it was only 380,000 years 
old and the transitional phase of cosmic re-ionization at redshifts z ~ 15-6. This 
epoch is commonly called the "dark ages" since it preceded - or, more precisely, 
included the beginning of - the formation of the first structures, which were either 
individual stars (or groups of stars) of the Population III, or an entire generation 
of stars within pre-galactic gas clouds. The radiation of these stellar objects lit 
up the sky and contributed significantly to the re-ionization of the intergalactic 
gas, a process that was basically completed by a redshift of z « 6 at an age of 
~ 800 million years. 

The first synthesis of the elements, mainly nucleons, occurred in a very short 
time interval approximately three minutes after the initial expansion of the uni- 
verse, but it took about four hundred thousand years before conditions were 
suitable for matter and radiation to decouple, with the former recombining into 
neutral atoms. However, recombination lagged behind expansion and at a red- 
shift of z ~ 800 a minute fraction of electrons and protons froze out at a level of 
one part in ~ 10 4 . The rapidly varying level of ionization, along with the ever 
decreasing energy of the photons of the cosmic microwave background (CMB) 
radiation, effectively promoted the formation of the first molecular species, ions 
and neutrals, that survived the processes of photoionization/photodissociation 
and dissociative recombination. This was, in fact, the "dawn of chemistry". 

Since Big Bang nucleosynthesis produced hydrogen and helium (along with 
their isotopes), with trace abundances of Li Be B and no solid particles, chemical 
reactions in the low density medium of the expanding universe only occurred in 
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the gas phase. Therefore, the first molecules and molecular ions were quite sim- 
ple and limited to few dominant species, such as H2, HD, HeH + , and LiH (Lepp 
&l Shull 1984). The search for these species prior to structure formation, either 
directly in the pregalactic medium or through their effects on the spectrum of 
the CMB has proved as yet unsuccessful in spite of several efforts (de Bernardis 
et al. 1993; Persson et al. 2010). Owing to its potential role in primordial chem- 
istry and in the cooling of the gas, LiH has also been searched for towards high 
redshift QSOs. Recent observations by Friedel, Kemball & Fields (2011) with 
the CARMA array of the J = — 1 transition of 7 LiH and 6 LiH have resulted 
in a 3cr-level detection of 7 LiH only in the source B0218+357 (z a bs = 0.685), 
confirming an earlier claim by Combes & Wiklind (1998). Similarly, a tentative 
detection of the J = 1 — rotational line of HeH + has been made towards one 
of the highest redshift quasars at z = 6.4 where other molecular lines had been 
previously detected (Zinchenko, Dubrovich &: Henkel 2011). 

Another obstacle in probing the physical conditions of the gas during the dark 
ages is that measurements towards high redshift galaxies and the diffuse gas have 
largely failed to identify clouds of primordial composition. On the contrary, there 
appears to be a systematic contamination by heavy elements at levels of Z ~ 10~ 3 
of the solar abundance (e.g., Wolfe, Gawiser & Prochaska 2005) which appears 
to extend up to z ~ 5 and which implies enrichment from the ejecta of earlier 
generations of stars (Rafelsky et al. 2012). However, recent studies have identified 
gas clouds without any trace of heavy elements or with metallicities consistent 
with enrichment by primordial stars (Fumagalli, O'Meara & Prochaska 2011; 
Cooke, Pettini & Steidel 2011). Interestingly, the pristine clouds LLS 1134a and 
LLS0956B with Z < 10" 4 Z Q at z « 3 identified by Fumagalli, O'Meara & 
Prochaska (2011) are located in regions of density significantly higher than the 
background, but unpolluted by heavy elements, indicating that mixing of the 
metals in the intergalactic medium is both inhomogeneous and inefficient. 

Finally, although a genuine zero-metal star has so far escaped detection, the 
search for such objects proceeds unabated due to the unique information that they 
carry about the earliest phases of galactic formation and evolution (e.g., Yong et 
al. 2012). The most exciting discovery of the last few years is the low- mass star of 
the galactic halo SDSS J102915+172927 with a metallicity 4.5 x 10~ 5 times that of 
the Sun and without enrichment of C, N, O (Caffau et al. 2011,2012). Although 
a number of stars of extremely low metallicity had been previously detected, 
they all showed significant overabundance of the heavy elements, a signature of 
pollution by a previous generation of massive stars (Frebel et al. 2008; Christlieb 
et al. 2002; Norris et al. 2007). From the analysis of their sample, Caffau et 
al. (2011) conclude that SDSS J102915 is not an isolated object and that a 
significant number of similar, or of even lower metallicity, low-mass stars will be 
discovered in the near future. These findings will definitely clarify the issue on 
the existence of a transition in the mass scale of the stars formed in the primordial 
universe. While a transition from massive to low-mass fragments based on the 
cooling properties of atomic lines of Ol, Cn, Sin and Fen and on the existence 
of a critical value of the C/H and O/H abundance has been suggested for some 
time (Bromm & Loeb 2003, Frebel, Johnson & Bromm 2007), the discovery of 
SDSS J102915 with metal abundance below the critical value appears to favor the 
view of a transition based on the cooling properties of dust formed in the ejecta 
of core-collapse SNe of the first generation (e.g., Schneider et al. 2012; Klessen, 
Glover & Clark 2012). 
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The exciting developments in the search for primordial molecules, clouds and 
stars outlined above make a review of the critical chemical processes that took 
place in the post-recombination universe timely. Until recently, large uncertain- 
ties existed not only on the rates of some key chemical reactions, but also (and, 
to some extent, more significantly) on the cosmological model and its fundamen- 
tal parameters, and on the yields from Big Bang nucleosynthesis (BBN). The 
advent of precision cosmology and a concordance model, along with significant 
improvements in standard BBN calculations offer accurate initial and boundary 
conditions for the study of the chemistry of the primordial gas, so that much of 
the remaining uncertainty resides in the poor knowledge of some basic reactions. 
Therefore, a coherent picture of the chemistry of the main constituents (hydrogen, 
deuterium, helium, and lithium) can be provided and, hopefully, tested. 

In this review, we will cover the post-recombination epoch in the unperturbed 
universe prior to the development of the first structures, whether individual stars 
or galaxies. Once these objects form, feedback effects and chemical enrichment 
become extremely important and need to be taken into account for a consis- 
tent treatment of the interstellar /intergalactic medium. Thus, the chemistry will 
be very different from that described here and an appropriate discussion would 
require a dedicated review. Similarly, the problem of the fragmentation and col- 
lapse of primordial clouds and the characteristic properties of the first stars will 
only be addressed in the light of the chemical properties of the gas. This field 
of research has made considerable progress in recent years thanks to the devel- 
opment of highly sophisticated numerical simulations and excellent summaries 
can be found in the literature (e.g., Bromm et al. 2009; Bromm &: Yoshida 2011; 
Loeb 2010; Umemura & Omukai 2013). 

The organization of this review is as follows: we sketch in Section [2] the early 
motivations to study the formation of molecular hydrogen and other molecules 
in the cosmological context. A good recipe for the chemical evolution of the 
primordial gas requires the specification of the main ingredients (a cosmological 
model, the abundance of the light elements from nucleosynthesis, the ionization 
history through recombination) and this will be provided in Section |3j Section [4] 
will present the results of the calculations of the predicted abundances of the 
main molecules and molecular ions, along with a discussion of the most critical 
reaction rates. The effects of chemistry on the cooling of the primordial gas will 
be highlighted in Section [5] Once formed, molecules as well as atoms can interact 
with the CMB and possibly leave some imprints both as spectral distortions and 
spatial anisotropics and Section [6] will address this important aspect. Finally, a 
summary of the major issues (Section[7| and some future perspectives (Section[8]) 
on both theoretical and observational studies of primordial chemistry will close 
the review. 

2 EARLY DEVELOPMENTS 

The investigation of the role of molecular hydrogen as an important coolant 
for the formation of the first structures in the early universe began in the late 
1960s. Saslaw & Zipoy (1967) were the first to point out the important role of 
H2 molecules for the thermal and dynamical evolution of pre-galactic gas clouds 
in the post-recombination era. Unlike previous attempts that had ignored the 
possible presence of H2 molecules due to the lack of dust grains and the slow rate 
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of radiative association and three body reactions, Saslaw & Zipoy considered the 
possibility of radiative association and charge transfer reactions such a^j] 

H + H + ^H+ + /u/ (1) 

and 

H+ + H^H 2 + H+. (2) 

The rate constant for reaction ([!]) was taken from a calculation by Bates (1951), 
while the cross section for reaction ^ was estimated to have a relatively high 
value of ~ 1CP 15 cm 2 , typical of charge exchange reactions above threshold. Thus, 
an ion can be converted to H2 as soon as it is formed and the H^" concen- 
tration remains very low. In addition, owing to the low electron abundance, 
dissociative recombination of turned out to be negligible compared with re- 
action ([2]) and H2 molecules could then be formed at a significant rate. Using the 
radiative cooling function calculated by Takayanagi & Nishimura (1960), Saslaw 
& Zipoy (1967) showed that, under reasonable assumptions for the initial condi- 
tions of the collapsing clouds, their evolution would depart from pure adiabatic 
collapse at densities above ~ 10 4 cm -3 due to enhanced H2 formation and cooling. 
Cloud contraction could then proceed at high densities before the onset of optical 
depth effects and H2 collisional dissociation, thus allowing the formation of stars 
of unspecified characteristics. Saslaw k, Zipoy suggested that the resulting object 
might resemble a globular cluster. 

Peebles & Dicke (1968) also considered the origin of globular clusters prior to 
the formation of galaxies. They noticed the correspondence between the char- 
acteristic mass and size of the first generation of bound systems, given by the 
Jeans mass and length (~ 10 5 Mq and ~ 5 pc, respectively), and that of globular 
clusters. However, Peebles & Dicke relied on another channel for efficient H2 
formation and cooling provided by the presence of negative hydrogen ions. They 
referred to McDowell (1961) who calculated the H2 abundance in the interstellar 
medium, following an earlier private communication by Dalgarno (19580 In this 
scheme, H2 formation proceeds through the reactions 

H + e H + hu (3) 

and 

FT + H -> H 2 + e. (4) 

The rate coefficient for ^ was calculated from the cross section given by Chan- 
drasekhar (1958), while that for reaction Q was based on laboratory measure- 
ments by Schmeltekopf, Fehsenfeld & Ferguson (1967). Once again, since the rate 
for Q is much more rapid than for ([3]), the formation of a negative hydrogen is 
promptly followed by the conversion to a hydrogen molecule. In addition, this 
route is faster than that through formation (by a factor of ~100 at 1000 K), 



1 They acknowledge P. Solomon, private communication 

2 Pagel (1959) was the first to point out the H~ channel for H2 formation in the solar photo- 
sphere, acknowledging A. Dalgarno for the suggestion. Later, McDowell (1961) calculated the 
H2 fractional abundance in Hi regions, quoting Pagel's calculation for the Sun and Dalgarno's 
estimate for the rate coefficient. Interestingly, McDowell mentions the fact that McCrea (1960) 
had pointed out that "consideration of the problem of star formation in HI regions strongly 
suggests that the hydrogen is present in molecular form". McDowell did not consider 3-body 
reactions because of the low densities, while McCrea & McNally (1960) suggested that the 
mechanism for H2 formation could be that of surface reactions on dust grains. 
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so that the latter was neglected in the calculations. The efficient formation of 
H2, up to a fraction of ~ 1CT 3 by number, and the enhanced cooling found by- 
Peebles Sz Dicke (1968) led to fragmentation of a small fraction of the cloud mass, 
and to collapse into stars mostly of high mass. 

The initial suggestions by Saslaw Sz Zipoy (1967) and Peebles & Dicke (1968) 
paved the route to more refined calculations of the role of H2 formation and cool- 
ing in collapsing and fragmenting primordial clouds, mainly with the goal of con- 
straining the typical masses of the forming objects in pregalactic gas clouds (e.g., 
Hirasawa, Aizu & Taketani 1969; Takeda, Sato & Matsuda 1969; Yoneyama 1970; 
Hutchins 1976; Silk 1977; Carlberg 1981). The general result was that even a rel- 
atively minor fraction of H2 contributes significant cooling via rotational and vi- 
brational transitions. However, owing to the lack of dipole moment, H2 molecules 
have large radiative lifetimes, thus are poor radiators. Nevertheless, these initial 
calculations showed that during collapse the gas temperature could be kept from 
climbing above ~ 1000 K over a large density range. As a result, the Jeans mass 
dropped by many decades and the minimum mass varied between ~ 20 M and 
~ 250 Mq with a maximum fractional abundance of H2 of about 10~ 3 before the 
onset of collisional dissociation. 

The inclusion in the chemical network of three body reactions 

H + H + H— >H 2 + H (5) 

and 

H + H + H 2 ^H 2 + H 2 . (6) 

modified this picture considerably (Palla, Salpeter & Stahler 1983). It was found 
that, over a wide range of initial conditions, virtually all the gas can be converted 
to molecular form by densities of the order of 10 12 cm -3 . As a result of significant 
cooling from the molecules, the temperature rise during collapse is slowed and 
the Jeans mass eventually drops below 0.1 Mq, thus allowing the formation of 
low-mass stars even in clouds of primordial composition. 

Much of the subsequent effort in the field of primordial chemistry was dedicated 
to elucidating the atomic and molecular processes, other than those related to 
H2, that are important in the early universe after the recombination epoch and 
prior to first structure formation and reionization. Major studies on this subject 
include: Lepp & Shull (1984), Black (1990), Puy et al. (1993), Standi, Lepp & 
Dalgarno (1996, 1998), Galli & Palla (1998, 2000, 2002), Puy & Signore (2002, 
2007), Lepp, Stancil & Dalgarno (2002), Glover & Abel (2008), Schleicher et 
al. (2008), Glover & Savin (2009), Vonlanthen et al. (2009), and Gay et al. (2011). 

3 THE INGREDIENTS OF PRIMORDIAL CHEMISTRY 
3.1 Cosmological Model 

Over the last few decades, the so-called A-Cold Dark Matter (ACDM) model 
has emerged as the standard description of cosmology. Accordingly, the uni- 
verse is spatially flat, homogeneous, and isotropic on large scales with a scale 
invariant spectrum of primordial seed perturbations composed of ordinary mat- 
ter, non-baryonic dark matter and dark energy. The accurate measurements of 
the peaks in the power spectrum of the CMB, along with baryon acoustic oscil- 
lations, the precise estimate of the Hubble constant, of light element abundances 
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Table 1: Cosmological parameters 
Parameter Value 



H 


100 h km s _i Mpc" 1 


h 


0.704 




3141 




10.6 


To 


2.725 K 




0.228 


n h 


0.0455 








^m/ (1 -2eq) 


n A 


0.727 


n K 




fu 


0.924 


fue 


0.076 


/d 


2.38 x 10~ 5 


fu 


4.04 x 10~ 10 



in distant quasars, and the amplitude of mass fluctuations inferred from clusters 
and gravitational lensing have enabled a deeper and more precise understanding 
of cosmology (e.g., Spergel et al. 2003; Komatsu et al. 2011). The essential pa- 



rameters (primary and inferred) of such a precision model are listed in Table 3.1 



All the values have been taken from the WMAP seven-year mean (cf. Table 1 
of Komatsu et al. 2011) with the exception of the current temperature of the 
CMB from Fixsen (2009), and the element abundances from Iocco et al. (2009). 
In the Table, z eq represents the redshift of equal energy density between matter 
and radiation, while z- lon is the redshift of reionization on the assumption that 
the universe was re-ionized instantaneously. The fractional abundances of H, He, 
D and Li relative to the total number of baryons are given by /(H), /(He), etc.. 
These cosmological parameters will be adopted in the rest of the discussion and 
will provide the initial values for the chemical evolutionary models that will be 
described in the next Section. It is expected that the Planck satellite by ESA 
with higher resolution, sensitivity, and frequency coverage than predecessors will 
provide an update on the current cosmological parameter estimates as early as 
2013. 



3.2 Standard Big Bang Nucleosynthesis 

According to ACDM cosmology, the formation of the first elements occurs entirely 
during the radiation-dominated epoch at a time t ~ (T/MeV) -2 s, where T is 
the plasma temperature. This energy scale corresponds to a time scale of a 
few minutes, temperature of ~ 10 9 K, and redshift of z ~ 10 8 -10 9 . The hot 
and low-density plasma was primarily composed of free nucleons and electrons 
that initiated the production of light elements (e.g., Fields 2011). Owing to 
the low density, most of the reactions occur by neutron and proton capture, 
whereas three body reactions are effectively suppressed, as well as reactions with 
nuclei heavier than ZiZj > 6. Thus, at temperatures T ~ 100 keV, deuterium 
(with binding energy 2.2 MeV) can be synthesized via p{n,^)D reactions, since 
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Table 2: Predicted Abundances from standard BBN 



Light Element PArthENoPE THALYS 





0.2467 


0.2476 


D/H (xl0~ 5 ) 


2.56 


2.59 


3 He/H (xl(T 5 ) 


1.02 


1.04 


7 Li/H (xlO- 10 ) 


4.60 


5.24 


6 Li/H (xlO- 14 ) 


1.11 


1.23 



blackbody photons can no longer equilibrate the reverse reaction. The abundance 
of deuterium rises sharply, as almost all available neutrons are locked into D- 
nuclei. The actual abundance is set by the value of the baryon-to-photon ratio, 
parametrized by the quantity 

rj = ^ = 2.74 x 10" 8 Q h h 2 (7) 
where fib = Phi Per , Pb is the baryon density, p CT is the critical density defined by 

and h is the present-day value of the Hubble constant measured in units of 
100 km s _1 Mpc -1 . In the following, we will adopt 7710 = 10~ 10 r/. 

The WMAP determination of r/10 = 6. 16 ±0.15 implies that the standard BBN 
has no free parameter left and can yield precise predictions of the abundances of 
D, 3 He, He, and 7 Li (Serpico, Esposito k, Iocco 2004; Iocco et al. 2009; Coc et 
al. 2012). These abundances can then be used as the initial conditions for the 
modeling of the chemical evolution of the expanding universe. Table [2] reports 
the values of BBN calculations obtained with two numerical codes: PArthENoPE 
(Iocco et al. 2009; Pisanti 2012, priv. comm.), and the nuclear reaction code 
THALYS (Coc et al. 2012). Although the codes can solve the synthesis of heavy 
elements (including CNO, 9 Be, and B), in the following we only consider the 
production of light elements since their abundance is sufficient for the chemistry 
calculations^] The predictions of the detailed models of BBN can be compared 
with the observations of light elements at high z. The literature on this subject 
is quite extensive and here we briefly summarize the current status (for a review, 
see, e.g., Steigman 2007). 

3.2.1 Deuterium. As for deuterium, there is consistency between the in- 
ferred primordial abundance deduced from the angular power spectrum of the 
CMB, and the predictions of BBN. Assuming that the measured D-abundance in 
low-metallicity environments at redshifts z > 2-3 corresponds to the true D/H 
ratio produced by the BBN (see, e.g., Kirkman et al. 2003), observations of 
Dl and Hi absorption lines in damped Ly-a (DLA) systems towards QSOs have 
provided the best determinations of deuterium. The mean value derived from 8 



3 However, the accurate assessment of the predicted abundances of CNO are of particular 
relevance for the evolution of Population III stars whose properties are significantly affected for 
values higher than ~ 10 -10 in mass fraction (Cassisi & Castellani 1993), or higher than 10 -12 for 
low-mass primordial stars (Ekstrom et al. 2010). This is still at least three orders of magnitude 
above the yields computed by Coc et al. (2012) and Iocco et al. (2009). 
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QSO absorbers is (Fumagalli, O'Meara & Prochaska 2011) 

log (D/H) p = (2.78 ±0.21) x 1(T 5 . (9) 

This implies a value of Q h h 2 (BBN)= 0.0213 ± 0.0012, in agreement within 
the errors with the value of ft^h 2 = 0.0222 ± 0.0004 deduced from the angu- 
lar power spectrum of the CMB (Keisler, Reichardt & Aird 2011). However, 
the authors note that the scatter between different measurements is larger than 
the published errors. The preferred explanation of the dispersion is an under- 
estimate of the systematic errors affecting D/H measurements (e.g., Pettini et 
al. 2008|5 Recen %> Pettini & Cooke (2012) have obtained an accurate measure- 
ment towards the metal-poor (Z ~ 1O~ 2 Z0), DLA at redshift z = 3.05 in the 
z cm ~ 3.03 QSO SDSS J1419+0829 with (D/H) p = (2.535 ±0.05) x 10~ 5 which 
implies tt h h 2 (BBN)= 0.0213 ± 0.0012. Similarly, Noterdaeme et al. (2012) have 
derived a primordial abundance of (D/H) p = 2.8^g'g x 10~ 5 , consistent with the 
predictions of BBN using the constraint on Q^h 2 from WMAP7. 

3.2.2 Helium. The situation for helium is less favorabl^B Since 4 He is 
produced in stars, the determination of its primordial value requires extrapo- 
lation from the measured abundances in metal-poor galaxies (extragalactic Hn 
regions), a procedure which is affected by systematic uncertainties (e.g., Olive 
& Skillman 2004, Izotov & Thuan 2010). Different analyses have yielded values 
that differ systematically from each other, with a weighted mean of 

Y p = 0.2566 ± 0.0023, (10) 

(Aver, Olive & Skillman 2010), that shows a slight discrepancy with the predicted 
BBN value reported in Table [2] We note that helium has the most accurately 
predicted abundances and it is quite insensitive to the precise value of the baryon 
density fib h 2 . In a recent re-analysis of the dataset of Izotov, Thuan & Stasin- 
ska (2007), Aver, Olive & Skillman (2012) using Markov Chain Monte Carlo tech- 
niques have obtained an extrapolated primordial value of Y p = 0.2534 ± 0.0083, 
in agreement with the WMAP7 result of 0.2487 ± 0.0002. 

3.2.3 Lithium. Finally, of all the light elements, lithium (in both isotopic 
versions, 7 Li and 6 Li) is most problematic, owing to the well established dis- 
crepancy between the predicted abundance of 7 Li and the value measured in 
metal-poor halo dwarf stars. In fact, the observations extended over a time span 
of about 30 years have confirmed the existence of a discrepancy between the mea- 
sured value of 7 Li/H= (1.58±0.31) x 10~ 10 - the so-called Spite plateau - and that 
inferred from WMAP measurements, larger by a factor of ~ 3 (Iocco et al. 2009). 
The realization of a dramatic scatter of the 7 Li abundance at very low metallici- 
ties, below [Fe/H] ~ —2.8, has complicated the interpretation of the evolution of 
lithium in low-mass stars even further (see Sbordone et al. 2010 and Bonifacio, 
Sbordone & Caffau 2012 for a recent analysis). A surprising, yet very encourag- 
ing, result has been reported by Howk et al. (2012) who have measured the 7 Li 
abundance in the interstellar medium of the Small Magellanic Cloud. Unlike the 

4 Alternative interpretations involving an early destruction of deuterium have been suggested 
(Olive et al. 2012). 

5 Since 3 He is both produced and destroyed in stars and since it can be observed only in the 
ISM of the Milky Way (e.g., Bania, Rood & Balser 2010), its abundance is affected by large 
uncertainties and cannot be used as a test for BBN predictions. 
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case for old metal-poor stars, the derived abundance of Li/H= (4.8±1.7) x 10~ 10 
is almost equal to the BBN predictions, thus providing a precious hint that the 
"lithium problem" posed by stellar measurements may not be representative of 
the primordial value. Since inferring the Li abundance in the ISM in not straight- 
forward and requires several assumptions, these initial results must be substan- 
tiated by further analysis. 

BBN predicts that 7 Li must be definitely more abundant than its isotopic ver- 
sion 6 Li by several orders of magnitude (cf. Table[2]) . High resolution observations 
of halo stars have been able to detect 6 Li with ratios as high as ~ 0.5 (Asplund 
et al. 2006; Garcia Perez et al. 2009). However, the observations remain con- 
troversial, as well as the evidence for a tendency of 6 Li to follow the same Spite 
plateau. Note that the observations towards the Small Magellanic Cloud give a 
limit to the isotopic ratio of 6 Li/ 7 Li < 0.28 (Howk et al. 2012). 

The standard BBN model rests on the assumption of homogeneity and isotropy 
of all constituent species, a hypothesis supported by the observed low amplitude 
of the CMB temperature fluctuations. However, alternative models of inhomoge- 
neous nucleosynthesis have been developed, mainly in order to provide a solution 
to the Li problem. In this case, the the baryon-to-photon and the neutron-to- 
proton ratios are allowed to vary on small scales allowing for a large variation 
of the element abundances (see, e.g., Iocco et al. 2009 for a general discussion). 
Here, it suffices to underline the potential impact that these models have on the 
yields of CNO and heavier elements as a distinctive diagnostic between standard 
and inhomogeneous models. Vonlanthen et al. (2009) have investigated the ef- 
fects of the higher CNO abundances resulting from non-standard BBN and their 
implications on molecule formation. 

3.3 Cosmological Recombination 

The importance of cosmological recombination for studies of primordial chem- 
istry and structure formation was already emphasized in the seminal works of 
Peebles (1968) and Peebles & Dicke (1968). According to Peebles (1968), "one 
also would like to know the residual ionization in the matter, because the electrons 
and ions mediate the formation of molecular hydrogen and so help determine the 
rate of energy radiation by clouds that formed at the Jeans limit". Figure 1 
shows the cosmological recombination of H, D, He and Li. Of all the nuclei of 
interest, Li 3+ is the first to recombine at redshift z ~ 14, 000, 



Li 3+ + e Li 2+ + hv, 



(11) 



then Li 2+ recombines at z ~ 8600 




(12) 



followed by the recombination of He 2+ at z ~ 6000, 




(13) 



and of He+at z ps 2500, 



He + + e ->• He + hv, 



(14) 



By z ~ 2000 all these reactions have run to completion, owing to the high density 
of the matter and the large availability of electrons as long as H is still ionized. 
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Figure 1: Cosmological recombination of H, D, He, and Li. Note the small resid- 
ual fraction of H + (and He + ), while D + is completely removed at low redshifts 
by charge exchange with H. Unlike all other elements, lithium remains largely 
ionized even at the lowest redshifts. 

Finally, at z ~ 1300 hydrogen recombines, along with deuterium: 

H+ + e^H + /izA (15) 

Notice that the cosmological recombination of Li, 

Li+ + e -> Li + hu, (16) 

essentially never occurs, owing to the photoionization of Li due to ultraviolet 
photons produced by the recombination of H and He (Switzer & Hirata 2005) 
and the five orders of magnitude drop of electron abundance following H recom- 
bination. Notice also that the residual fraction of D + , unlike those of H + and 
He + , does not reach a freeze-out because D + is efficiently removed at low redshift 
by the near-resonant charge-transfer reaction 

D+ + H^D + H+, (17) 

which is exothermic by 3.7 meV (43 K). 

3.3.1 Recombination of Hydrogen. The basic physics of cosmological 
H recombination was fully understood in the late 1960s (Peebles 1968; Zel'dovich 
et al. 1968; see Peebles 1993). Direct recombination to the n = 1 ground state 
of H or He produces photons which can easily ionize neighboring neutral atoms, 
because the primordial gas has a large optical depth (~ 10 7 ) to ionizing photons. 
Therefore, cosmological H and He recombination proceeds by capture of electrons 
into excited n > 2 states followed by decay to the ground state. However, elec- 
trons in excited states cannot cascade instantaneously down to the ground state, 



12 



D. Galli, F. Palla 



because of the large reservoir of photons in the CMB with energies sufficient to 
produce reionizations out of these states (E > 3.4 eV). In fact, the population 
of these excited states is, with good approximation, in thermal equilibrium with 
the CMB. In addition, transitions to the ground state produce Lyman photons 
that can excite other neutral atoms allowing them to be easily re-ionized. Unlike 
the present-day ISM, where Lyman photons are eventually absorbed by dust, the 
only way they can escape reabsorption in the conditions of the early universe is 
through two processes of low efficiency: by redshifting out of the line because of 
the cosmological expansion, or by two-photon forbidden emission from the n = 2 
level. The "bottleneck" caused by slow decay from the n = 2 to the n = 1 
level, combined with the density drop due to the expansion, prevents H recom- 
bination from being completed and results in a residual fraction of electrons and 
protons. On the basis of a simple two-level approximation (n = 1, n > 2) plus 
continuum, Peebles (1968) and Zel'dovich et al. (1968) found a residual ioniza- 
tion 2 x 10 -5 -2 x 10~ 4 depending on the present-day value of the mass density, 
in agreement with previous estimates of Novikov & Zel'dovich (1967). Many 
refinements and additional effects were later added to the original formulation 
(e.g. Matsuda, Sato & Takeda 1971; Jones & Wise 1985; Krolik 1990; Sasaki & 
Takahara 1993) without however changing the basic concepts summarized above. 
Still in 2000, the major uncertainty on the residual electron fraction was due to 
the poor knowledge of the cosmological parameters fib and h. 

As a result of the exquisite accuracy of experiments such as WMAP, several de- 
tailed multilevel calculations, supplemented by computer codes, are now available 
to determine with precision better than 0.1% the evolution of the electron frac- 
tion down to z » 500: RECFASlQ (Seager, Sasselov & Scott 1999,2000; updated 
by Wong, Moss & Scott 2008 and Scott & Moss 2009); RICC[](Fendt et al. 2009); 
COSMOREC0 (Switzer & Hirata 2008ab; Ali-Haimoud & Hirata 2010; Chluba, 
Vasil & Dursi 2010; Grin & Hirata 2010; Rubiho-Martm et al. 2010; Chluba & 
Thomas 2011); HyRecQ (Ali-Haimoud & Hirata 2011); ATLANT0 (Kholupenko 
et al. 2011). 

3.3.2 Recombination of Helium. Cosmological He recombination was 
first computed by Matsuda, Sato & Takeda (1969). An updated and detailed 
review of the processes involved can be found in Switzer & Hirata (2008a, b,c) 
and Chluba, Fung & Switzer (2012). In currently available computer codes like 
RECFAST, H and He recombinations are computed simultaneously, with He + 
and He 2+ treated as multilevel atoms. The recombination of He 2+ presents no 
n = 2 "bottleneck" because of the high two-photon decay rate. Consequently, 
the abundance of He 2+ is well approximated by the Saha formula. This is not 
true for He + recombination, which follows essentially a hydrogenic-like case B 
recombination and is found to be slow essentially for the same reasons that H 
recombination is slow (see, e.g., Seager, Sasselov & Scott 2000). 

3.3.3 Recombination of Lithium. The calculation of Li recombination 
has not been straightforward. It was first computed by Palla, Galli & Silk (1995), 
and later refined by Standi, Lepp & Dalgarno (1996), Bougleux h Galli (1997) 

6 http://www. astro. ubc.ca/people/scott/recfast. html 

7 ht tp :// cosmos . astro . uiuc.edu/rico 

8 http://www.cita.utoronto.ca/~jchluba/ Science- Jens/Recombination/CosmoRec. html 
9 http://tapir. caltech.edu/~yacine/hyrec/hyrec. html 

10 http://www.ioffe.ru/astro/QC/CMBR/atlant/atlant.html 
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and Galli & Palla (1998). Denning the recombination redshift z re c,Li as the red- 
shift when Li/ (Li + Li + ) = 0.1, Stancil, Lepp & Dalgarno (1996) and Galli 
& Palla (1998) found Zrec,Li ~ 400—450, and a residual Li ionization fraction 
Li + /Li ~ 0.5-0.7 as z — > 0. However, it was later realized by Switzer & Hi- 
rata (2005) that the recombination history of Li is strongly affected by UV pho- 
tons generated by the recombination of H and He at z < 1000 via Ly-a and 



two- photon decay (see Section 3.3.1). These nonthermal photons are sufficient 
to reduce the neutral Li fraction by about 3 orders of magnitude at z = 100. 
The resulting evolution with redshift in Figure 1 shows that Li never recombines 



completely and Li + is still ~ 98% of the total Li at z = 10 (cf. Table 4.1 ). As dis 



cussed in Section 6.1, the recombination history of Li has important cosmological 
consequences. 

4 CHEMICAL EVOLUTION AFTER RECOMBINATION 

The full network of chemical reactions for the early universe has become increas- 
ingly complex over the last decade, from the ~ 90 reactions of Galli & Palla (1998) 
to the ~ 250 reactions of Gay et al. (2011). The number of species formed from 
nuclei up to Li is limited to ~ 30, but for state-resolved chemistry each state must 
be considered as a separate species. Thus, for example, H^~ and H2 require 19 and 
15 coupled equations, respectively, one for each vibrational state of the funda- 
mental electronic level (Coppola et al. 2011). The network has been extended to 
molecular species produced by elements heavier than Li, like F (Puy et al. 2007) 
and C, N, O (Vonlanthen et al. 2009), but their abundances are extremely small 
even if non-standard BBN models are adopted. Ignoring the very small change 
in particle number due to chemical reactions, the total baryon density evolves 
with redshift as 

n b = ^(l + .) 3 (18) 

where p CT is the critical density defined by equation ^ , mn is the mass of the 
hydrogen atom, and \x = 4/(4 — 3Yp) the mean atomic weight of the primordial 
gas. Numerically, ~ 2.2 x 10~ 7 (1 + z) 3 cm~ 3 . 

Once formed, molecules can be photodissociated by the CMB photons with 
rates that depend on its spectrum, a blackbody characterized by a radiation 



temperature T r = Tq(1 + z), with To = 2.725 K (see Table 3.1). In addition, 
distortions of the CMB produced by photons emitted during the recombination 
of H and He cannot be neglected, as they affect significantly, and in some cases 
dramatically, the rates of several photodestruction processes, such as H~ pho- 



todetachment (see Section 4.1) and Li recombination (see Section 3.3.3). 

For applications to the homogeneous universe, the chemical reaction network 
must be supplemented by a routine that provides the ionization fraction as func- 
tion of redshift for the adopted cosmological model, and an equation for the 
evolution of gas temperature T g , 

AT R 2 

= ~^TgJ^ + 3fc nfa ~~ A)compton + (r — A) mo i + Tchem]- (19) 

In this equation, the first term represents the adiabatic cooling of the gas due to 
the expansion of the universe, R being the scale factor; the second and the third 
the net transfer of energy (heating T minus cooling A) from the radiation to the 
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gas via Compton scattering of CMB photons on electrons, and via excitation and 
de-excitation of molecular transitions, respectively; the fourth term represents 
chemical heating of the gas (see e.g. Glassgold, Galli & Padovani 2012). The 
ionization fraction is sufficient to maintain T g ~ T r up to z ~ 300 via Compton 
heating of the gas. Because of the small molecular abundances, the term (T — 
A)moi is very small and can raise the gas temperature by only a few degrees 
(Khersonskii 1986; Galli & Palla 1989; Puy et al. 1993). The chemical heating 
term r c h em is also negligible (Puy et al. 1993). Consequently, at redshift z < 
100, the gas cools almost adiabatically, with T g ~ 0.02(1 + z) 2 K. Additional 
sources of heating for the gas that could affect its chemical evolution are dark 
matter particles decays and annihilations (Ripamonti, Mapelli & Ferrara 2007) 
and ambipolar diffusion of a cosmological magnetic field (Schleicher et al. 2009). 
Finally, the chemical network is completed by the equation for the redshift, 

dt __ 1 

Tz ~ ~{i + z)H{ z y ( ' 

where H{z) is the Hubble parameter 

H(z) = H [n T (i + zf + n m (i + zf + n K (i + zf + 0\] 1/2 . (21) 

The meaning of these quantities has been discussed in Section |3.1| and their 
numerical values are listed in Table |3~T1 



4.1 Overview of Global Trends 

The evolution of the main molecular species using the cosmological parameters 
and initial element abundance described before is displayed in Figure 2 as frac- 
tional abundances relative to the total number of baryons. It can be conveniently 
described following the five main periods during which some common process can 
be identified that dominates the chemical kinetics of the gas. As a function of 
decreasing redshift, the sequence is as follows: 

I. 2000 < z < 800. At z ~ 1000 the ionization fraction is still ~ 10%, while at 
z = 800 it drops to ~ 1%. While hydrogen is still mostly ionized, the availabil- 
ity of protons promotes the initial steep rise of all molecular species. The first 
reaction to occur is the radiative association of H + and He 

H+ + He -)• HeH + + hu, (22) 

quickly followed by the sequence of radiative association of H and H + (reaction 
[I]) that leads to H2 formation via charge exchange (reaction [2]) . The formation 
of H2" via reaction Q is more efficient than the conversion of HeH + to H^ by 
reaction with hydrogen atoms. In this phase, the main destruction process is 
through energetic photons from the CMB. As shown in Figure 2, the formation 
of H2, HD, H^, and H2D + proceeds quite rapidly and then levels off at a redshift 
of z ~ 1000. However, the HeH + , H^, and HD + abundances display a turnover 
followed by a decrease to a minimum at z ~ 800 due to recombination with 
electrons. 

II. 800 < z < 300. As the ionization fraction freezes out, all the molecular 
species start forming at a moderate pace, sustained by the decreased impor- 
tance of photodissociation processes by the CMB. The major limiting factor to 
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Figure 2: Fractional abundances (relative to the total number of baryons) of the 
main molecules and ions formed in the early universe as function of redshift. The 
vertical dotted lines indicate the boundaries of the five main evolutionary phases 
described in the text. 



the increase of H2 molecules is the importance of deviations from local thermo- 
dynamic equilibrium (LTE) in the H^ level population and the highly efficient 
photodissociation by CMB photons (Hirata & Padmanabhan 2006; Coppola et 
al. 2011). Thus, the behavior of HD, H^, and H2D + follows exactly that of the 
H2 molecules. In principle, the kinetics of HeH + formation should also be fol- 
lowed including the effects of non-LTE analysis of the level populations, but this 
calculation has not been done yet. One would then expect that the steep rise 
shown in Figure 2 could just reflect the assumption of a LTE level population and 
that the actual abundance should be considerably reduced. At redshifts z < 400, 
LiH begins to form by radiative association, but its growth is limited by efficient 
photodestruction. 

III. 300 < z < 60. The abundance of most species rises and reaches a peak 
value at z ~ 60. As the energy density of the CMB photons keeps decreasing, 
forms more effectively through radiative association of H and H + , but the charge 
exchange reaction [2] limits its abundance. At redshifts below ~ 100, the formation 
of H ions starts becoming important as the main route for H2 formation. H ions 
(binding energy 0.754 eV) are no longer effectively photodetached by the CMB 
and promote the formation of H2 molecules by radiative attachment (reaction [4]) . 
However, as first computed by Hirata & Padmanabhan (2006), the presence of 
a significant tail of non-thermal photons of the CMB produced by hydrogen and 
helium recombination slows down the rise of H - abundance and, hence, of H2 
and H3". 

IV. 60 < z < 20. The abundance of H2 and HD reaches freeze-out, while 
that of H steadily drops due to its conversion into H2 molecules. Freeze-out 
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Table 3: Elemental and molecular abundances relative to hydrogen at z = 10 



H 2 /H 
6.3 x 10~ 7 


H"/H 
1.8 x 10~ 13 


H+/H 
9.2 x 10~ 15 


H+/H 
8.0 x 10~ 17 


H 3 -/H 
5.5 x 10~ 33 


HD/H 


HD+/H 


D+/H 


H 2 D+/H 




4.2 x 10~ 10 


1.2 x 10~ 18 


6.1 x 10~ 19 


1.6 x 10" 20 




HeH+ /H 
1.7 x 10~ 14 


He+/H 
9.8 x 10~ 23 


He+/H 
9.3 x 10~ 36 






Li+/H 
4.3 x 10~ 10 


Li/H 
8.0 x 10~ 12 


LiH/H 

9.0 x 10~ 20 


LiH+/H 
4.6 x 10" 20 


Lr/H LiHe+/H 
1.7 x 10" 22 3.0 x 10" 23 



is achieved when the formation time becomes longer then the Hubble expansion 
time. At redshift z < 40 the marked increase in the HeH + and HD + abundance 
results from the exponential drop in the reaction 

HeH + + H -> He + H+ (23) 

below T g = 100 K (Bovino et al. 2011c). Some of the HeH + participates in the 
enhanced formation of HD + via 

HeH + + D -»• HD+ + H, (24) 

which drives further production of H^ and HD + . Below z ~ 40, LiH + is readily 
formed due to the decrease of the intensity of the CMB, and its abundance reaches 
that of LiH. 

V. z < 20. As first structure formation begins to be important, the steady de- 
crease of some species is accounted for by the relevant role played by dissociative 
recombinations. Interestingly, only H 2 and HD have reached their steady state 
abundances, whereas the rest remain out of equilibrium. This fact introduces 
some arbitrariness in the initial conditions for the calculations of the dynam- 
ical evolution of the density perturbations, even though the minute fractional 
abundances of most chemical species do not affect the thermal properties of the 
primordial gas. 

The chemical evolution of each element will now be described in some detail. 
4.2 The Chemistry of Hydrogen 

4.2.1 H 2 formation THROUGH excited states. During the recombina- 
tion era, the formation of H 2 molecules can proceed either via radiative associa- 
tion of excited hydrogen atoms (Latter &; Black 1991) 

H(n = 2) + H(n=l)^H 2 (A 1 S+) + /i^ (25) 

or via Raman association (the inverse Solomon process; Dalgarno & van der Loo 
2006; Alizadeh & Hirata 2011) 

H + H + /«/^H 2 (5 1 S+,C 1 n u )^H 2 (X 1 S+) + ^, (26) 

where A 1 !!^ is the ground state of H 2 , and -B 1 E+ and C l Ii u are the excited 
electronic states with energies ~10 eV. The potential curves of the and 
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C 1 H U with allowed transitions to the X 1 ^^ ground state of H2 correlate with 
atomic hydrogen in n = 2 (both 2s and 2p) and n = 1 levels, respectively. 

In the first process, as pointed out by Latter &i Black (1991), the population of 
the n = 2 excited level is maintained over the LTE value as a result of resonant 
scattering by Ly-a photons (Peebles 1968; Zel'dovich et al. 1968), thus allowing 
reaction (25) to occur rapidly at a typical rate of few 10 -14 cm 3 s _1 over a 
temperature interval of 50-30000 K. Latter & Black showed a pronounced increase 
of H2 molecules due to this reaction in the redshift interval between 2500 and 
1500. 

Recombination of H2 by Raman association is facilitated owing to the abun- 
dance of ultraviolet photons, both blackbody CMB photons (Dalgarno & van der 
Loo 2006) and spectral distortion photons (Alizadeh & Hirata 2011). The result- 
ing excited H2 molecule reemits the photon and can decay either to the bound 
state X 1 S+ or to the continuum of the X level leading to dissociation. Dalgarno 
& van der Loo (2006) computed the H2 production rates via Raman association 
including all relevant rovibrational resonances and found them to be more domi- 
nant than the radiative association channel of reaction (25 ). In these calculations, 
most of the contribution to H2 formation occurs at a redshift z ~ 1600. 

On the other hand, Alizadeh & Hirata (2011) have demonstrate the importance 
of the distortion photons of the CMB using a code that included the absorption 
and emission of Lyman-series photons in the Hi lines. This contribution starts to 
be dominant at redshifts z < 1400 when the photon density deviates significantly 
from the Planck spectrum. The overall effect is a net reduction of the H2 abun- 
dance due to efficiency of the inverse of reaction (25). In the redshift interval 
1300-900, the H2 fractional abundance remains well below the value of 10 -16 , 
several orders of magnitude lower than the value computed neglecting this pro- 
cess. At lower redshifts, however, the process becomes rapidly negligible as the 
spectral distortion is redshifted at energies lower than the Lyman- Werner bands. 
Interestingly, Alizadeh & Hirata (2011) find that the H2 abundance is too low 
to affect the recombination process or to introduce secondary anisotropics in the 
CMB: in principle, absorption or emission of H2 molecules of photons above the 
Lyman- Werner bands could add or remove photons that would otherwise excite 
a hydrogen atom and alter the recombination history. 

4.2.2 State-resolved formation of H 2 and H^. At redshifts lower 
than z ~ 800, the abundance of H2 molecules is determined by the balance be- 
tween charge transfer with H^~ ions, reaction Q and the reverse of reaction |l|. 
Other minor routes involve reactions of HD molecules with either H or H + . In 
spite of its simplicity, the kinetics of the main reactions for H2 is quite complex 
due to the presence of a large number of quantum states and the evaluation of 
the exact vibrational distribution function is indeed necessary (see, e.g., Lepp, 
Stancil & Dalgarno 2002). This aspect has proved particularly important in the 
case of the photodissociation rates of Hg" ions, the main destruction route in com- 
petition with conversion into H2 molecules, which have usually been computed 
either assuming LTE level populations of the rovibrational levels or with all the 
H^ ions in the ground state. Detailed calculations of the fully level-resolved 
populations of H^ have been performed by Hirata & Padmanabhan (2006) and 
Coppola et al. (2011; see also Longo et al. 2011). The resulting behavior is shown 
in Figure 3a. The main conclusion of these studies is that full thermalization of 
the vibrational levels is never achieved since the vibrational relaxation processes 
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are too slow to balance the strong vibrational selectivity of formation rates and 
radiative processes affect both the vibrational distribution function and the over- 
all fractional abundances. As a consequence of the non-LTE level population, 
the formation of H2 molecules is greatly reduced mainly because the H^" ions 
are photodissociated before they can either decay to the ground state or undergo 
charge transfer to become H2 molecules. 

4.2.3 H 2 formation via H~ ions. The last step in the formation of H 2 
molecules is via associative detachment, reaction Q. Recently, the cross section 
for this simple reaction has been accurately measured in experiments with merged 
beams apparatus at collision energies of cosmological interest (Kreckel et al. 2010; 
Miller et al. 2011). Once formed, H~ ions are readily destroyed by photons of 
the CMB owing to the rather low photodetachment threshold energy (0.754 eV). 
The major improvement on the precise magnitude of this process has been the 
inclusion of the spectral distortions in the CMB (Hirata & Padmanabhan 2006). 
The rate coefficient for photodetachment can be written as 

fcnonth = nnc / r v o v — (27) 

where r v is the number of recombination photons per H atom and per logarithmic 
frequency range, u v is the photodetachment cross section and hv§ = 0.754 eV. 
The nonthermal photons become the dominant process at redshifts below z ~ 130 
when the photon energy is ~ 1 eV, corresponding to the peak of the photode- 
tachment cross section. The effect on the abundance is clearly visible in 
Figure 3b where the sharp peak at z ~ 100 (dotted line) is replaced by a shal- 
low increase and a maximum at z ~ 60. This behavior accounts for the evolution 
of H 2 molecules, also shown in Figure 3b. Unlike the case without non-thermal 
photons, the fractional abundance of H2 reaches the freeze-out value only at 
z < 60, when the nonthermal photons have been redshifted below the threshold. 
The freeze-out value of H 2 molecules at z = 10 is 6 x 10~ 7 . 

4.2.4 Minor species: H~, and H^. As shown in Figure 3b and Ta- 



ble 4.1, the fractional abundance of these species at z = 10 is very small, varying 



between ~ 10~ 13 for H and ~ 10~ 16 for H^. The reason why H~ ions do not 
reach equilibrium at the lowest redshifts is that they are constantly removed by 
the radiative association with H atoms to form H2 molecules. Similarly, the 
ion displays a steady decrease from its peak abundance of ~ 2.5 x 10~ 16 at z « 50 
due to efficient destruction by dissociative recombination with electrons. Accord- 
ing to the experimental results of McCall et al. (2004) and Kreckel et al. (2005), 
the dissociative reaction rate is in general agreement with the theoretical predic- 
tions and the branching ratios tend to favor the channel H+H+H over the H+H2 
channel. The global evolution of follows quite closely that of H2 as the main 
formation process is via radiative association with H + at all redshifts. Due to its 
potential importance as a coolant of the collapsing gas, we shall return to in 



the discussion on the rates (Section 4.6). 



As for H2 , the behavior of is strongly affected by the inclusion of the non- 
LTE level population in the calculation of the photodissociation rate, which is 
the main destruction process of this ion. The latter ceases to be important at 
redshifts below z ~ 300 and that explains the flattening of the abundance shown 
in Figure 3a. At lower redshifts, formation by radiative association is balanced 
by charge exchange reaction into H 2 , accounting for a nearly constant fractional 
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abundance of ~ 10 . Finally, at redshifts smaller than z ~ 60, the formation of 
H^~ is enhanced by the HeH + channel via the reaction with whose rate drops 
significantly at temperatures below T g ~ 100 K (Bovino et al. 2011c), as will be 
discussed in Section 14.41 

We also mention the possibility of the formation of Hg , the negatively charged 
counterpart of H^", which is predicted to be stable by about 0.013 eV (Starck 
k, Meyer 1993). Recently, Ayouz et al. (2011) have developed the theory of its 
formation in the ISM by radiative association of H2 and H - and derived the 
relevant cross sections. Using the reaction rate by Ayouz et al. (2011), we have 
verified that the resulting abundance is neg ligible (< 10 -30 ) at all redshifts. 



4.3 The Chemistry of Deuterium 

Deuterium in the early universe is important for a variety of reasons. First, ow- 
ing to its larger mass, deuterated hydrogen has a spacing of the rotational levels 
smaller than that of hydrogen molecules. Thus, cooling through HD dipole radi- 
ation can effectively drive the primordial gas to substantially low temperatures, 
provided a sufficient abundance is obtained. Thus, the chemistry of deuterium 
has been computed by numerous groups following the initial study of Lepp & 
Shull (1984), such as Dalgarno & Lepp (1987), Latter (1989), Puy et al. (1993), 
Palla, Galli & Silk (1995), Standi, Lepp & Dalgarno (1998), Galli & Palla (1998, 
2002), Lepp, Standi & Dalgarno (2002), and Glover & Abel (2008). Recently, 
Gay et al. (2011) have presented a thorough analysis of the chemistry of highly 
deuterated species, including D2, D^ , D2H + , and Djj~. The formation of HD 
follows two main routes involving a deuteron exchange with H2 

D + H 2 -> HD + H, (28) 

and 

D + + H 2 ^HD + H+. (29) 

The thermal rate constant for reaction ( |28| ) has been measured in numerous exper- 
iments over a wide range of temperatures (e.g., Mitchell & LeRoy 1973; Michael 
k, Fisher 1990); theoretical calculations (e.g., Mielke et al. 1994; Charutz, Last & 
Baer 1997) show very good agreement with each other and with the experimen- 



tal data. Reaction ( 29 ) represents the major source of HD in diffuse interstellar 
clouds (Dalgarno, Weisheit & Black 1973). Its rate coefficient is almost constant 
and close to the Langevin value 2.1 x 10 -9 cm 3 s _1 (e.g., Gerlich 1982). Finally, 
at z < 10, HD formation takes place via charge exchange (Karpas, Anicich & 
Huntress 1979) 

HD + + H — >• HD + H + . (30) 



The destruction of HD occurs mainly via the reverse reactions of (28) and (29) 

HD + H -> D + H 2 (31) 

and 

HD + H+^D+ + H 2 . (32) 
The rate coefficient for the first reaction has been computed by Shavitt (1959) 



using a semi-empirical H3 energy surface. As for reaction (32), owing to its 
endothermicity by 0.0398 eV (462 K), the removal of HD at low temperatures 
is reduced by a factor exp(— 462 K/T g ), leading to significant fractionation. At 
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redshift (1+z) redshift (1+z) 

Figure 3: (a) Hydrogen chemistry. The solid curves show the case with state- 
resolved formation of H2 and Hjj", while the dashed lines correspond to the evo- 
lution obtained by ignoring this effect. The dotted lines display the evolution 
of H2 and H~ without the contribution of the non-thermal photons of the CMB 
to the photodissociation of H~. (Adapted from Coppola et al. 2011). (b) Deu- 
terium chemistry. The solid and dashed curves have the same meaning as in 
the caseAstron. Astrophys. of the hydrogen chemistry shown in Figure 3a. 
(Adapted from Galli &: Palla 2002). (c) Helium chemistry. Evolution of the main 
He-bearing molecular ions, HeH + and He^. (Adapted from Bovino et al. 2011c). 
(d) Lithium chemistry. Evolution of Li-bearing molecules and ions with (solid 
curves) and without (dotted curves) the effects of non-thermal photons on Li 
recombination. In the latter case, Li never recombines, molecular abundances 
decrease significantly at all redshifts, with the exception of Li - , and LiH remains 
the most abundant species. (Adapted from Bovino, Tacconi & Gianturco 2011b 
and Bovino et al. 2011d). 
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freeze-out, the ratio HD/H2 ~ 7 x 10~ 4 , with an enhancement factor of ~ 25 over 
the initial D/H abundance. 

The relative abundance of D and D + , is governed by the charge exchange 
reactions 

H+ + D^H + D+, (33) 

and 

H + D+ -> H + + D. (34) 



Since reaction (33) has a threshold of 43 K, the rate coefficient for reaction 



(34) is usually obtained by multiplying that of reaction (33) by exp(43 K/T g ). 
Savin (2002) has computed accurate rate coefficients for both reactions, im- 
proving on the early estimates by Watson (1976) and Watson, Christensen & 
Deissler (1978) that have been widely used in studies of deuterium chemistry. 

Figure 3b highlights the close similarity of the evolution with redshift of the 
HD abundance with that of H2 and the dramatic effect of including the state 
resolved levels in the calculations. Generally, the chemistry of HD is dominated 



by the ion-neutral reactions ([29]) and (32) in a gas of low density (e.g., prior to 



the formation of the first structures), whereas the neutral- neutral reactions (28) 



and (31) become more important in conditions of high density and temperature 
(as during cloud collapse or in shocked gas). 

Figure 3b also displays the evolution of the main deuterated cations HD + and 
H2D + . The case of H2D + is of interest due to the suggestion of Dubrovich (1993) 
and Dubrovich & Lipovka (1995) that its abundance could reach significant val- 
ues (up to H2D + /H2 = 10 -5 ), in which case specific rovibrational transitions 
could leave detectable imprints on the spectrum of the CMB. In principle, such 
an abundance could result from the complete conversion of into H2D + fol- 
lowing deuterium fractionation at low temperatures (T g <C 50 K). However, the 
actual abundance of H3" in the primordial gas is quite low due to fast dissociative 
recombination. Although Dubrovich & Lipovka (1995) claim observable features 
in the CMB for abundances of H2D + as low as 10 -8 , Figure 3b shows that this 
value is never achieved in the post-recombination universe where H2D + reaches 
a peak of only ~ 2 x 10 -19 at z ~ 80. On the other hand, the isotopic ratio is 
enhanced with respect to the primordial D/H abundance by a factor of ~ 8. As in 
the case of H^j", the main removal of H2D + occurs via dissociative recombination. 

Finally, the chemistry of HD + is directly related to that of Hjj" since its abun- 
dance is determined by the equilibrium between the forward and reverse reaction 

H+ + D -> HD+ + H. (35) 

HD + would maintain the freeze-out value of ~ 5 x 10~ 19 , reached at a redshift 
of ~ 250 if not for the enhancement due to the production from HeH + via 

HeH+ + D -> HD+ + He. (36) 

Comparison of Figure 3a and Figure 3b visualizes the similar behavior of HD + 
and H2", as well as H2D + and Hg~. 

4.4 The Chemistry of Helium 

The helium hydride ion (HeH + ) and the molecular helium ion (HeJ) contend 
for the role of the first molecular species formed in the universe. Searches for 
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HeH + in space, either in gaseous nebulae (Moorhead et al. 1988, Liu et al. 1997, 
Dinerstein & Geballe 2001) or in high-redshift absorbers (Zinchenko, Dubrovich 
&: Henkel 2011), have only given tentative detections. 

In the early universe, HeH + and He^" are formed at z ~ 7000 by radiative 
association of He with H + and He + , respectively. The rate of the former process is 
well determined theoretically (Saha, Datta h, Barua 1978; Flower & Roueff 1979; 
Roberge & Dalgarno 1982; Basu & Barua 1984; Kimura et al. 1993; Dumitriu & 
Saenz 2009; Sodoga et al. 2009) and the absolute photodissociation cross section 
has been measured by the free electron laser FLASH (Pedersen et al. 2007), 
although at high photon energy (38.7 eV). Some studies have also addressed the 
enhancement of the photodissociation rate produced by CMB photons (Jufek, 
Spirko & Kraemer 1995; Zygelman, Stancil &: Dalgarno 1998). The rate of the 
He^ formation has a temperature dependence similar to that of the radiative 
association of H and H + , but is about one order of magnitude smaller (Stancil, 
Babb & Dalgarno 1993). 

The history of helium chemistry is displayed in Figure 3c. Since He + re- 
combines rapidly at z ~ 2500, the abundance of He^ reaches a maximum of 
~ 10~ 23 around this redshfit, and then is readily destroyed by photodissociation 
and dissociative recombination. In contrast, the incomplete recombination of H + 
keeps the abundance of HeH + steadily rising to ~ 10~ 13 at z » 20. The mam 
reactions responsible for the bumpy behavior of the evolution of HeH + and that 
limit its final abundance are photodissociation for z > 300, collisions with H for 
20 < z < 300 (Linder, Janev & Botero 1995, Bovino et al. 2011c), and dissocia- 



tive recombination for z < 20 (Guberman 1994). As mentioned in Section 4.2 
destruction of HeH + by collisions with H atoms is an important source of H^ and 
accounts for its enhancement at the smallest redshifts. Finally, at redshifts below 
z « 800, photodissociation processes become inefficient, and the abundance of 
He^~ rises again by about 5 orders of magnitudes, but its growth is limited by 
efficient destruction by collisions with H for z < 400. 

In spite of the potential importance of HeH + molecules for efficiently scattering 
the CMB photons (see Section [6]) , an accurate calculation of photodissociation 
rates is still lacking. Miyake, Gay &: Stancil (2011) have obtained cross sections 
for all rotational transitions from the vibrational levels v" = 0-11 of the 1 XT, + 
state for temperatures between 500 and 12,000 K, including both X X and 
A X transitions. These cross sections can then be used to directly obtain the 
photodissociation rate due to the CMB radiation. Miyake et al. present some ini- 
tial estimate for different radiation temperatures and find that at redshifts > 1000 
the A <— X transitions (which are usually neglected) dominate the photodisso- 
ciation, while at lower redshifts the X X transition become more important. 
Thus, one would expect that the initial rapid rise of the HeH + abundance shown 
in Figure 3c will be further reduced by inclusion of this photodestruction channel. 
Future work is needed in order to clarify the magnitude of this effect. 

4.5 The Chemistry of Lithium 

The chemistry of lithium in the early universe has been reviewed by Dalgarno, 
Kirby & Stancil (1996), Stancil, Lepp & Dalgarno (1996), Bougleux & Galli (1997), 
Galli & Palla (1998), Puy & Signore (2002) and Bodo, Gianturco & Marti- 



nazzo (2003). After the recombination of Li 3+ and Li 2+ (see Section 3.3), the 



main remaining species Li and Li + react with the residual electrons, H, D, and 
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He atoms and ions to form Li - , LiH, LiH + , LiD, LiD + , LiHe + and Li - . The 
chemical network of lithium is complex (for updated chemical reaction rates see 
Bovino, Stoecklin & Gianturco 2010; Bovino et al. 2010; Bovino, Tacconi & 
Gianturco 2011b; Bovino et al. 201 Id; Bovino et al. 2012), but the molecular 
abundances are very small. In addition, as these species are easily photodissoci- 
ated, their formation in significant amounts is delayed to relatively low redshifts, 
z < 200. To complicate matters, the evolution is also extremely sensitive to the 
effects of nonthermal photons from the recombination of H and He (cf. Switzer 
h Hirata 2005). The latter can easily maintain Li ionized down to the lowest 
redshifts, whereas in their absence lithium would recombine at the same level of 
Li + . This behavior is displayed in Figure 3d where the dramatic effect of the 
inclusion of non-thermal photons can be appreciated by comparing the evolution 
in the two cases. In addition to the neutral and singly ionized species, some Li~ 
ions can also form at redshifts below z ~ 100. Their evolution follows closely 
that of H - , both anions being removed at low temperatures mainly by reactions 
of mutual neutralization with H + , with a rate typically proportional to T g . 
As a result, the abundance of Li - at z = 10 is ~ 2 x 10~ 22 . 

The first Li-bearing molecule to form is LiH mainly via radiative association, 

Li + H -> LiH + hi/, (37) 

with a contribution for 20 < z < 80 from Li~ via the reaction 

LP + H -»• LiH + e. (38) 

Note, however, that the rate for the latter reaction is only estimated. Radiative 
association takes place first from the Li( 2 P) state, then from the Li( 2 5) state. 
Since the 2 P level lies about 1.85 eV above the 2 S level, the formation of LiH 
molecules via electronically excited levels is important only at z > 800, when the 
temperature is greater than some thousand degrees K. Early estimates for the 
radiative association rate from the 2 S state resulted in relatively large values, 
and therefore the abundance of LiH in the early universe was expected to be 
a substantial fraction of the Li abundance (Lepp Sz Shull 1984). However, the 
situation was clarified by refined quantum mechanical calculations that reduced 
significantly the magnitude of the cross sections (Dalgarno, Kirby & Stancil 1996; 
Gianturco & Gori Giorgi 1996). As for destruction, LiH is easily photodissociated 
by CMB photons for z > 300, whereas at lower redshifts it is mainly destroyed by 
collisions with H, a process for which quantal calculations are available (Bovino, 
Wernli & Gianturco 2009). 

Owing to the small binding energy (0.14 eV), the main production of LiH + 
via radiative association of Li with H + is delayed until redshifts below z ~ 40, 
although a minute contribution can result at z > 700 from reactions of Li + ions 
with H. The reaction rates for both channels have been computed with fully quan- 
tal methods (Dalgarno, Kirby & Stancil 1996; Gianturco & Gori Giorgi 1997). 
As in the case of LiH, photodissociation is the dominant destruction mechanism 
of the molecular ion, with an uncertain contribution from collisions with H and a 
contribution from dissociative recombination (Curik & Greene 2007,2008) at very 
low redshifts (z < 20). Interestingly, at z = 10 the abundance ratio of LiH/LiH + 
is ~ 2, either including or ignoring the non-thermal photons of the CMB. 

Finally, LiHe + is formed by spontaneous and stimulated radiative association 
of Li + with He (Bovino, Tacconi & Gianturco 2011a, 2012). At z > 100, LiHe + is 
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destroyed by photo-dissociation, while at lower redshift by collisions with H and 
by dissociative recombination (Bovino et al. 2012). However, the abundance of 
LiHe + at z = 10 remains exceedingly small, ~4x 10~ 23 . 

4.6 Critical Rates 

As stressed in the Introduction, the evolution described so far depends more on 
the uncertainties of the chemical model than on the initial conditions provided 
by the cosmological models. The abundances of the various species are only as 
accurate as the rates given as input. The literature is rich in studies that critically 
examine the uncertainties of the reaction rates and their effects on the chemistry 
(and cooling) of the primordial gas. Most of the important reactions involved in 
the chemistry of the major species have well determined rate coefficients. Below, 
we give a brief overview of the rates for the most relevant reactions or that still 
present significant differences either from theory or from experiments. Extensive 
and updated compilations of reaction rates can be found in Glover & Abel (2008), 
Schleicher et al. (2008), Glover & Savin (2009), Bovino et al. (201 Id), and Gay 
et al. (2011). 

4.6.1 Associative detachment of H~. Surprisingly enough, in spite of 
many years of research, there has been considerable uncertainty on both the mag- 
nitude and temperature dependence of the rate of this important reaction that 
affected all the calculations of the chemical and thermodynamical evolution of 
Population III stars (e.g., Glover, Savin & Jappsen 2006). The series of merged 
beams measurements carried out by Bruhns et al. (2010), Kreckel et al. (2010), 
Miller et al. (2011) and the ion trap study by Gerlich et al. (2012) have given 
results in agreement with the theoretical predictions (e.g., Sakimoto 1989, Ci'zek, 
Horacek k, Domcke 1998). Figure 4 illustrates the good match (to within ±25%) 
between experimental data of Gerlich et al. (squares), Bruhns et al. (solid line) 
and theory (dashed curve) over a wide temperature interval and highlights the 
departure from the Langevin value that had been previously adopted in chemical 
models. The impact of the newly computed rates on the collapse and fragmen- 
tation of primordial clouds has been discussed by Kreckel et al. (2010) who have 
shown that the uncertainty on the predicted Jeans mass of the fragments is re- 
duced from a factor of 20 to only a factor of 2 (see Savin et al. 2012; Savin 2013). 

4.6.2 Mutual neutralization of H and H+. Although this process is 
irrelevant in the chemistry of the expanding universe, it becomes important when 
the gas gets reionized and recombines out of equilibrium during the gravitational 
collapse phase (see, e.g., Glover, Savin k, Jappsen 2006). The older experiment 
by Moseley, Aberth & Peterson (1970) reports cross sections that are consistently 
higher at energies above 3 eV than the most recent ones by Szucs et al. (1984) 
and Peart &: Hayton (1992). The latter are in very good agreement with the the- 
oretical calculations (see, e.g., Stenrup, Larson & Elander 2009). Experiments at 
energies below ~ 3 eV are required in order to remove the remaining uncertainty 
(Urbain et al. 2010). 

4.6.3 Radiative association of This reaction which dominates the 
kinetics is important at all temperatures and redshifts. The rate coefficient 

calculated by Stancil, Babb & Dalgarno (1993) is in excellent agreement with 
the older results of Ramaker & Peek (1976). Below ~ 10 K, the rate is basically 
constant at the value of ~ 2 x 10~ 20 cm 3 s _1 which accounts for the flat portion 
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Figure 4: Rate coefficients (in cm 3 s _1 ) for some critical reactions. Upper left 
panel: Filled squares: data from Gerlich et al. (2012) with la statistical un- 
certainty. Black line: results from the merged-beam experiments of Bruhns et 
al. (2010) and Miller et al. (2011). Red curve: calculation by Cfzek et al. (1998). 
Triangles: measurements at room temperature by Schmeltekopf, Fehsenfeld 
&z Ferguson (1967), Fehsenfeld, Howard & Ferguson (1973), and Martinez et 
al. (2009, slightly shifted for better readability). Short-dashed line: Langevin 
value. Upper right panel: Blue and red curves: fits by Galli &; Palla (1998) and 
Stancil, Lepp & Dalgarno (1998) of the rate by Schneider et al. (1994) for the 
v = 0, J = 0-20 levels (triangles). Black curve: LTE rate computed by Coppola 
et al. (2011) from the state-resolved cross sections of Takagi (2002). Lower left 
panel: Triangle: experimental result at 300 K (Karpas, Anicich & Huntress 1979). 
Red and blue curves: reaction rate obtained by Stancil, Lepp & Dalgarno (1998) 
and Schleicher et al. (2008) from the experimental cross sections of Linder, Janev 
& Botero (1995). Black curve: result of the quantum calculations by Bovino 
et al. (2011c). Lower right panel: Li( 2 S'): data from Dalgarno, Kirby & Stan- 
cil (1996, cyan), Gianturco & Gori Giorgi (1996, blue), Bacchus-Montabonel h 
Talbi (1999, red), Bennett et al. (2003, 2008, green). Li( 2 P): data from Gi- 
anturco & Gori Giorgi (1997, blue) for the B X IL -> X l Y> + (filled dots) and the 
A l Y, + — > (empty dots), respectively. 
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of the curve of the H^~ abundance at low redshifts (see Figure 3a). We note 
that the fitting formulae given by Stancil, Lepp & Dalgarno (1998) and Gay et 
al. (2011) become inaccurate below T g ~ 10 K, possibly affecting the resulting 
H^" abundance at the lowest redshifts. 

4.6.4 Photodissociation of H^. This is the main destruction process 
for which determines the history of H2 molecules down to redshifts z ~ 
300. State-resolved cross sections have been calculated by Dunn (1968), Stancil, 
Babb & Dalgarno (1993), Lebedev, Presnyakov & Sobel'Man (2000), and Hirata 
&: Padmanabhan (2006). However, only the early results of Dunn (1968) are 
available in tabular form and have been adopted in the work by Coppola et 
al. (2011) on which the discussion of Section 4.2 is based. 

4.6.5 Charge exchange of H 2 with H + . This endoergic reaction has 
a threshold of 1.83 eV. The most complete theoretical work on this chemical 
reaction is that of Krstic (2002). The thermal rate coefficient has been computed 
by Savin et al. (2004) for R 2 (v = 0, J = 0). Glover & Abel (2008) discuss the 
large uncertainties on this rate that exist in the literature, quoting as extreme 
values the prescription by Shapiro &: Kang (1987) based on detailed balance from 
the reverse reaction, and the much lower rate of Abel et al. (1997), based on cross 
sections by Janev, Langer & Evans (1987). However, the latter is valid only at 
temperatures higher than ~ 10 4 K, and should not be used at lower temperatures 
where it exponentially drops to unrealistic values. At temperatures above ~ 
10 3 K, the rate coefficient by Savin et al. (2004) is in excellent agreement with 
the fit by Galli & Palla (1998) based on cross sections of Holliday, Muckerman 
&: Friedman (1971). At lower temperatures, the exponential drop due to the 
energy threshold makes the uncertainty in the rate coefficient of no consequence 
for primordial chemistry and collapse calculations. 

4.6.6 Charge exchange of H2 WITH H. Also for this exoergic reaction, 
state-resolved cross sections have recently been computed by Krstic (2002) for 
collision energies in the range 0.6-9.5 eV, but need to be extrapolated to lower 
energies for applications to primordial chemistry. The extrapolation is usually 
made in such a way as to reproduce the value of the thermal rate of this reaction 
measured by Karpas, Anicich & Huntress (1979) at T g = 300 K (Coppola et 
al. 2011). Clearly, specific calculations at low energy (and additional experimental 
measurements) are needed. 

4.6.7 Dissociative recombination of H^. A wealth of theoretical and 
experimental data is available for the dissociative recombination of H^" and its 
isotopologues (see, e.g., Schneider k, Suzor-Weiner 2002). The level of exper- 
imental accuracy made possible by the advent of storage-ring experiments has 
encouraged a continuous refinement of the computational techniques and the in- 
clusion of previously neglected physical processes. Using multi-quantum defect 
theory, Takagi (2002) has calculated state-resolved dissociative recombination 
cross sections that include the rotational structure and all excited dissociative 
states. A relevant characteristic is the presence of many resonances at low ener- 
gies, and the sensitivity of the cross section to the vibrational quantum number 
v, especially at low energies. As a result, the reaction rate computed by Cop- 
pola et al. (2012) in LTE has a non-monotonic behavior significantly different 
from the usual power-law approximations commonly adopted in the literature 
(see Figure 4). 
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4.6.8 Three-body H2 formation. In the low density gas of the expand- 
ing universe, three-body reactions cannot take place. However, in the collapsing 
clouds of the primordial halos where the first generation of stars is thought to 
occur, these reactions effectively transform the atomic gas into fully molecular 
hydrogen at densities above ~ 10 9 cm -3 . In addition, since each reaction lib- 
erates 4.48 eV of energy, corresponding to the binding energy of H2, into the 
gas, three-body reactions represent an important energy source for the dynami- 
cal evolution of the clouds. The main reaction for H2 formation involves atomic 
hydrogen as the third body, but a number of other channels are also possible. 
Although there is consensus on the small value of the rate coefficient of the reac- 
tion, the numerical values quoted by different authors vary by up to two orders 
of magnitude at the low temperatures of interest for cosmological simulations, 
while there is substantial agreement (to within a factor of ~2) at temperatures 
greater than ~ 2000 K (see, e.g., Glover & Savin 2009). 

Current estimates of the three body reaction rates are based on laboratory 
measurements of the reverse reaction, collisional dissociation of H2 by hydrogen 
atoms, using the principle of detailed balance (e.g., Palla, Salpeter & Stahler 1983; 
Flower &: Harris 2007). Similar, although smaller, differences exist in the rate 
coefficients for reactions where the third body is an H2 molecule. The source of 
the discrepancy is attributed by Flower Sz Harris (2007) to different assumptions 
made for the equilibrium constant. However, including the nuclear degeneracy, as 
well as the electron spin degeneracy, in the calculation removes the discrepancy 
and reduces the rate to the values adopted by Palla, Salpeter & Stahler (1983) 
(F. Esposito, priv. comm.). An even lower value of the rate was proposed by Abel 
et al. (2002), based on the calculations of Orel (1987) at temperatures < 300 K 
and assuming an inverse temperature dependence at higher temperatures. From 
the theoretical viewpoint, Esposito & Capitelli (2009) have presented a unifying 
approach that includes the contributions from quasibound and unbound states 
by means of the classical orbiting resonance theory (e.g., Orel 1987) and from 
continuum states through directly calculated dissociation with detailed balance 
application. Their total recombination rate shows a dependence on the temper- 
ature to pressure ratio, but it is in fairly good agreement with the estimates of 
Abel, Bryan & Norman (2002) at low temperatures. 

Performing direct experiments of three-body association has proved extremely 
difficult in the laboratory due to the high densities and temperatures required. 
Savin (2013) gives an account of the challenges of such experiments and concludes 
that the situation is not too promising in the short term. Turk et al. (2011) have 
performed an analysis of the effects of the current uncertainty in the rate using 
3D simulations of the collapse of primordial clouds in the high density regime. 
The results clearly show how sensitive the outcome is to the precise value and 
temperature dependence on the rate, concluding that "the uncertainty represents 
a major limitation on our ability to accurately simulate the formation of the first 
stars in the universe". Clearly, this is a field that urgently requires some basic 
technological or methodological breakthrough. 

4.6.9 Radiative association of Hj. Although the reaction 

H 2 + H + -> H+ + hv (39) 

is the main route for H^j" formation, the rate coefficient has an uncertainty of four 
orders of magnitude. The study of Gerlich & Horning (1992) based on an ion 
trap measurement and a classical trajectory analysis recommends the constant 
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value of 1 x 1CP 16 cm 3 s _1 . This value has been adopted in most chemical 
models (e.g., Schleicher et al. 2008; Glover & Savin 2009). However, Stancil, 
Lepp & Dalgarno (1998), on the basis of quantal calculation of other diatomic 
molecules, argue that the rate should be much smaller and indicate a value of 
1 x 10~ 20 cm 3 s -1 also adopted in the recent compilation of Gay et al. (2011). 
The evolution of Hjj~ is clearly affected by such a large uncertainty. If instead of 
the large rate, the value suggested by Stancil, Lepp & Dalgarno (1998) had been 
adopted in the chemical model discussed in Section 4.2 the resulting abundance 
of would decrease by three orders of magnitude with respect to that shown in 
Figure 3a. Thus, this reaction needs to be further investigated before the effective 
role of Hj|~ can be elucidated. 

4.6.10 Dissociative recombination Hg and H 2 D + . The ion storage 
ring measurements of McCall et al. (2004) and Kreckel et al. (2005) have settled 
the issue of the rate coefficient of this most relevant reaction. Previously, a large 
uncertainty existed between the results of merged beam experiments with typical 
values of ~ 10~ 7 cm 3 s _1 at room temperatures (e.g., Sundstrom, Mowat & 
Danared 1994) and of flowing afterglow experiments with values about ten times 
smaller (e.g., Smith & Spanel 1993). More recent experiments (e.g., Glosfk et 
al. 2008), as well as refined theoretical calculations have converged towards the 
results of McCall et al. (2004) with the remaining disagreement below a factor of 
2. In the evolution presented here, we have followed McCall et al. (2004) using 
a branching ratio of 0.66 for the recombination into H + H + H and of 0.34 into 
H2 + H + . As for H2D + , the experiments of Datz et al. (1995) and Larsson et 
al. (1996) still remain the main source for the determination of the rate coefficient 
with branching ratios of 0.73 for the channel H + H + D, 0.20 for HD + H, and 
0.07 for H 2 + D. 

4.6.11 Isotopic exchange OF HD. In spite of being the most important 
channel for the destruction of HD, only sparse laboratory data exist for reaction 
(31). In addition, all the experiments are limited to temperatures greater than 
~ 200 K (e.g., Galli &; Palla 2002), thus confirming the plea made long ago by 
Shavitt (1959) for more experimental work on such a basic reaction. 

4.6.12 Transfer reaction of HeH + . Reaction (23) is the main destruc- 
tion channel of HeH + in the redshift interval 300 < z < 10. Prior to Bovino et 
al. (2011c), it was not studied with sufficient accuracy from a theoretical point 
of view. The crossed-beam data for energies higher than 0.2 eV (Rutherford &; 
Vroom 1973) were fitted by Linder, Janev & Botero (1995) and adopted in the 
chemical models of Schleicher et al. (2008). Bovino et al. (2011c) present new ab 
initio quantum calculations for this process for collision energies from 10 -4 eV to 
1 eV. The reaction rate is displayed in Figure 4. The marked drop of the thermal 
rate below T g ~ 100 K is due to the non-Langevin behavior of the cross sections 
at energies below ~ 10 -2 eV. 

4.6.13 Radiative ASSOCIATION of LiH. The cross sections and rates for 
the radiative association from the Li( 2 P) state have been computed by Gianturco 
Sz Gori Giorgi (1997). As for the 2 S state, using semi-classical arguments, Lepp 
&z Shull (1984) estimated a radiative association rate ~ 10~ 17 cm 3 s" 1 , a value 
confirmed by Khersonskii & Lipovka (1993). Subsequent quantal calculations 
performed by two groups on the basis of independent sets of potential curves and 
wavefunctions (Dalgarno, Kirby &: Stancil 1996; Gianturco & Gori Giorgi 1996) 
provided a value of the rate smaller than the previous results by about three orders 
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of magnitude. Semiclassical calculations by Talbi & Bacchus-Montabonel (1998) 
and Bacchus-Montabonel & Talbi (1999) were found in agreement with the quan- 
tal results. The inclusion of quasi-bound levels with radiative widths exceeding 
their tunneling widths (Bennett et al. 2003) resulted in an increase of the rate 
by a factor of ~ 2. The rates of radiative association of LiH from the 2 S and 2 P 
states are shown in Figure 4 as function of temperature. 

5 COOLING OF THE PRIMORDIAL GAS 

Cooling and thermal balance of the primordial gas are key ingredients in the for- 
mation of the first structures, as they determine the characteristic fragmentation 
scale of the clouds and the magnitude of the accretion flow onto the growing pro- 
tostars (e.g., Abel et al 2002; Yoshida et al. 2006). While chemistry determines 
the abundances of the main molecular constituents of the gas, the microphysics 
of each individual component defines the cooling efficiency of the gas. In the case 
of the homogeneous expanding universe, the thermal balance of the gas is simply 
set by adiabatic expansion since molecules do not form in sufficient fraction to 
effectively couple the matter and radiation fields (see Section [4|. However, in 
the collapsing gas of the first minihalos or larger halos, the density and temper- 
ature increase considerably allowing the excitation of more molecular degrees of 
freedom. Thus, a detailed knowledge of the main cooling agents is needed to 
determine the thermal state of the gas. 

5.1 Molecular Cooling 

In order to evaluate the cooling function at each redshift, the radiative and colli- 
sional transition probabilities and the frequencies of the rovibrational transitions 
must be specified for each molecule. These data can be found in the literature 
not only for the most relevant molecules, like H2 and HD, but also for the other 
species, such as Hj", H^~, HD + , H2D + , HeH + , LiH, and LiH + . Then, for each 
molecule, one must solve the level populations from the rate equations. The 
studies by Glover & Abel (2008) and Glover & Savin (2009) provide detailed dis- 
cussions for each species. More recently, Coppola, Lodi k, Tennyson (2011) has 
computed updated radiative cooling functions valid in the high-density regime 
where LTE conditions prevail. One should notice that accurate collisional co- 
efficients are available only in limited cases, and for the minor species simple 
approximations are used in the evaluation of the cooling functions. We summa- 
rize in Figure 5 the individual cooling functions for the most relevant species in 
the low-density case (a) and in LTE conditions (6), respectively. 
5.1.1 H2 COOLING. In spite of its dominant role, H2 has two basic limita- 
tions: the excited ro- and vibrational transitions have small radiative transition 
probabilities (radiative lifetimes > 10 6 s) and collisional deexcitation becomes 
competitive with spontaneous transitions at rather low densities (n ~ 10 4 cm -3 ). 
Thus, as the density increase the cooling rate of H2 reaches the LTE value which 
depends only on the transition probabilities and the energy levels. Since the for- 
mer are small, the cooling rate is also small. The second effect is the onset of 
optical depth in the cores of the rovibrational lines that effectively reduces the 
cooling efficiency and requires the solution of the appropriate radiative transfer 
equations. Detailed calculations have shown that the cooling by H2 molecules 
is effectively suppressed at densities n > 10 10 cm~ 3 (e.g., Yoshida et al. 2006). 
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Figure 5: Cooling rate/function for H2, H^", H^j~, HD, and LiH in the low-density 
(a) and LTE (6) limits. In the low-density case, we plot the cooling rate of ortho- 
H2, para-H2 (dashed curves), and total H2 (solid curve) assuming an ortho/para 
ratio of 3:1. Sources for the individual low-density cooling rates are: Glover &; 
Abel (2008) for H 2 ; Suchkov & Shchekinov (1978) and Glover & Savin (2009) for 
Hj; Glover & Savin (2009) for H^; Lipovka, Nunez-Lopez & Avila-Reese (2005) 
for HD; Bougleux & Galli (1998) for LiH. Sources for the LTE cooling func- 
tions are: Coppola et al. (2012) for H 2 ; Glover k Abel (2008) for H^ ; Glover & 
Savin (2009) for H+; Coppola et al. (2011) for HD and LiH. 



Another important aspect is the H2 ortho/para ratio, since the equilibrium value 
of 3:1 is achieved only at relatively high temperature and densities. Indeed, in 
the expanding low-density universe, the H2 formation process can decrease the 
LTE value for abundances fn 2 < 10 -6 . Glover Sz Abel (2008) have shown that 
for H + /H> 10 and temperatures lower than several hundred K, the ortho/para 
ratio is controlled by collisions with protons, so that the cooling rate is mainly 
due to para-H2 (J = 2-0, AE/k = 510 K), although most of the gas contains 
ortho-H 2 (J = 3-1, AE/k = 845 K). 

In the low-density limit, the total cooling is the sum of individual contribu- 
tions due to the collisional partners of H2: H, H2, He, electrons, and protons. 
Glover & Abel (2008) provide a thorough discussion and empirical prescriptions 
of the various terms. Some recent improvements concern the H2-H2 collisions, 
discussed by Lee et al. (2008, see also Balakrishnan et al. 2011) who provide 
updated collisional coefficients based on the potential energy surface of Diep h 
Johnson (2000). 

5.1.2 HD cooling. Due to small but non zero dipole moment (8.3 x 
10~ 4 debyes), HD has radiative lifetimes that are a factor of ~ 100 shorter 
than those of H2, so that the cooling rate reaches its LTE limit only at den- 
sities n ~ 10 6 cm~ 3 . In addition, unlike H2 molecules, transitions with A J = 1 
are permitted, allowing cooling to occur through the J = 1-0 transition (corre- 
sponding to AE/k = 128 K). Finally, due to the exothermicity of reaction (29) by 
462 K, HD is highly fractionated at low temperatures, thus boosting the relevance 
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of this molecule as an important coolant for the gas at the lowest temperatures. 
The cooling function computed by Lipovka, Nunez-Lopez & Avila-Reese (2005) 
for HD-H collisions for 1 < n < 10 8 cm -3 and 100 < T„ < 10 4 K, still remains 
the main reference. In fact, it is in very good agreement with the results of 
Wrathmall, Gusdorf & Flower (2007) who have improved the interaction poten- 
tial and included a larger set of rovibrational levels. In the LTE limit, Coppola, 
Lodi & Tennyson (2011) find a slightly larger (factor of ~ 2) cooling rate for 
temperatures T g > 700 K due to the contribution of highly excited rovibrational 
levels. Collisions of HD with He and electrons have been studied by Nolte et 
al. (2012) and Yoon et al. (2010), respectively. However, to our knowledge, the 
corresponding cooling functions have yet to be evaluated. 

5.1.3 LiH and LiH + cooling. The very large dipole moment (5.89 debyes) 
of LiH makes this molecules an excellent radiator, and potentially important as a 
coolant of the primordial gas. However, its low abundance renders its contribution 
to the total cooling rate insignificant in the expanding universe while it may play 
a role in collapse calculations. The cooling function for collisions with H has 
been computed by Bougleux & Galli (1997, fitted by Galli & Palla 1998). These 
results are valid below a critical density of n CT ~ 10 12 cm -3 (Lepp & Shull 1984), 
and are still commonly used in the literature. At densities larger than n cr , the 
radiative cooling rate is provided by Coppola, Lodi & Tennyson (2011) who 
computed updated potential energy and dipole moment curves. As for LiH + , 
whose abundance is slightly smaller than that of LiH, only the LTE cooling 
rate is available (Coppola, Lodi & Tennyson 2011). Its value is a factor of 10- 
100 smaller than that of LiH, making LiH + a negligible coolant in cosmological 
conditions. 

5.1.4 Minor species. Cooling by H^ has been the subject of a detailed 
analysis by Glover & Savin (2009) and Miller et al. (2010). The former find 
that H^~ is the third important coolant after H2 and HD in dense (n ~ 10 7 - 
10 9 cm 3 ) primordial gas. In this density regime, the consequences of the l arge 



4.6 



uncertainty on the rate of radiative association of H3 discussed in Section 
are not as strong as those on the expanding medium, because other channels 
dominate the formation of this ion. As for H^ and HD + , Glover Sz Savin (2009) 
give cooling functions both in the low- and high-density regime. In the latter case, 
however, Coppola, Lodi &: Tennyson (2011) find a HD + cooling rate larger by 
orders of magnitude than that of Glover & Savin (2009) at temperatures < 500 K. 
Finally, the LTE cooling rate of HeH + and its isotopologues has been computed 
by Coppola, Lodi & Tennyson (2011) (see also Miyake & Standi 2007). Due to 
is large dipole moment (1.66 debyes, Pavanello et al. 2005), the LTE cooling rate 
is estimated to be about ten orders of magnitude grater than that of H2. No 
calculations are available for the low-density regime. 

5.2 Chemical and Thermal Evolution of Collapsing Primordial 
Clouds 

The well constrained initial conditions provided by the standard ACDM model 



(see Section 3.1 ) and the knowledge of the atomic and molecular processes in the 
primordial gas allow the investigation of the evolution of dark matter structure 
and the formation of the first objects. Current consensus is that the first genera- 
tion of stars started to take place at redshifts z ~ 20-30 in dark matter halos of 
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Figure 6: Thermal [upper) and chemical [lower) evolution of collapsing primor- 
dial gas clouds. The solid lines refer to the case of collapse in a predominantly 
neutral minihalo, while the dashed lines display the evolution of gas which is 
pre-ionized and cools from high temperatures. Lines of constant Jeans mass are 
shown as dotted lines in the upper panel. (Courtesy of K. Omukai). 



mass ~ 10 6 M.q (called minihalos) and preceded the formation of the first galaxies 
(e.g. Bromm 2012). In these minihalos, star formation could proceed thanks to 
the cooling provided by H2 molecules (Couchman & Rees 1986; Tegmark et al. 
1997) in a predominantly atomic gas and without the complications of physical 
processes (viz., turbulence, magnetic fields, etc.) that characterize the environ- 
ment of present-day molecular clouds (however, for a different view see, e.g., 
Schleicher et al. 2009; Clark et al. 2011; Turk et al. 2012). 

Along with the prevailing minihalos, the ACDM model predicts the presence of 
more rare, massive halos with mass ~ 10 8 Mq and characteristic temperature of 
~ 10 4 K, too hot for the survival of molecular hydrogen. In this case, it is believed 
that the gas can cool via atomic hydrogen (Oh & Haiman 2002), but fragmenta- 
tion is effectively suppressed with the possible formation of massive/supermassive 
black holes (Bromm &i Loeb 2003; Hosokawa et al. 2012). However, a situation 
in which molecular cooling becomes relevant is when the gas is reionized and re- 
combines out of equilibrium, as in expanding Hn regions and intergalactic shocks 
(e.g. Shapiro & Kang 1987). The resulting enhanced ionization fraction promotes 
the efficient formation of H2 via the H channel and HD molecules through D + 
via reaction (29). If the HD abundance becomes larger than a critical value, the 
gas can cool below the floor temperature of ~ 200 K produced by H2, and even 
approach the minimum value set by the CMB (Johnson & Bromm 2006; Yoshida, 
Omukai & Hernquist 2007; McGreer & Bryan 2008). 
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The thermal and chemical evolution of the collapsing primordial gas can be 
understood following the results of simplified 1-D calculations that capture the 
essential microphysics. Figure 6 illustrates the two cases discussed above: a 
predominantly neutral medium in a minihalo and a pre-ionized medium. 

5.2.1 Formation of first stars in minihalos. As shown in Figure 6 
(solid curves), due to the low initial fractional H2 abundance (/h 2 < 10 -6 ), 
the gas collapses almost adiabatically and H2 formation proceeds mainly via the 
H~ channel. At relatively low densities (n « 10 cm -3 ), the fractional abundance 
rises to fu 2 ~ 10~ 3 which is sufficient to greatly enhance the cooling efficiency via 
rovibrational transitions that effectively contrasts the gravitational compression 
of the gas and produces the marked drop in temperature shown in Figure 6 
up to densities n ~ 10 4 cm -3 . The minimum temperature of ~ 200 K sets a 
characteristic mass for fragmentation of 

/ n N-i/2 / T \ 3//2 

This mass scale is a robust result, rather insensitive to variations in the initial 
conditions. At higher densities, however, H2 molecules reach LTE conditions 
and line cooling saturates. The gas temperature rises slowly with density, while 
the fractional H2 and HD abundances remain basically constant at a value of 
/h 2 ~ 10 -3 and /hd ~ 10 -7 . The major change occurs at densities higher than 
~ 10 s cm -3 when three body reactions begin to convert the prevailing atomic 
hydrogen into H2 molecules, which partly offset the rise in temperature. The 
increased H2 abundance drives the accelerated formation of HD molecules. At 
densities higher than ~ 10 12 cm ~ 3 , the cooling is governed by H2 collision- 
induced emission and dissociation. Once H2 molecules are almost completely 
destroyed, the gas temperature increases adiabatically and pressure forces even- 
tually succeed in effectively halting gravitational collapse and a protostellar core 
is formed that accretes rapidly matter from the surrounding cloud. 

The final outcome of the collapse and the typical masses of the first stars are 
still subject to uncertainty. While for a certain time it seemed that the first 
stars were predominantly very massive with characteristic masses > 100 M 
(e.g., Abel et al. 2002; Bromm, Coppi & Larson 2002; Yoshida et al. 2006), 
more recent calculations including feedback effects have shown that the growth 
of protostars at the very high accretion rates (M > 10~ 3 M yr _1 ) expected due 
to the lack of coolants can be limited to < 50 M for a wide range of initial 
conditions (McKee & Tan 2008; Hosokawa et al. 2011). High-resolution 3-D 
simulations have also shown that multiple fragments can develop during collapse 
and circumstellar disk formation, thus leading to the formation of binary or small 
multiple systems or even small clusters of stars of typically low mass (Clark, 
Glover k Klessen 2008; Clark, Glover & Smith 2011; Turk, Abel & O'Shea 2009; 
Stacy, Greif & Bromm 2010). Considering the sensitivity of the accretion history 
on the thermal and dynamical state of the parent gas cloud and on the physical 
characteristics of the minihalos, it is likely that collapse and fragmentation of 
the primordial gas will produce a first generation of stars with a substantial 
distribution in mass. 

5.2.2 Formation of stars from pre-ionized gas. The thermo-chemical 
evolution of primordial gas (see Figure 6, dashed curves) differs significantly from 
that described above, especially in the initial phases of rapid cooling and enhanced 
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H2 formation that results from out-of-equilibrium recombination of electrons with 
protons. In these conditions, the higher ionization fraction promotes the rapid 
formation of H2 molecules that reach fractional abundances of about 5 x 10~ 3 at 
densities < 1 cm~ 3 , a factor of > 50 higher than in the neutral case. As a response, 
the gas drops to temperatures well below ~ 150 K, the threshold value for HD 
formation (Omukai et al. 2005). More HD molecules are readily formed, up to 
/hd ~ 10 -5 , while the gas cools to ~ 30 K at densities ~ 10 cm -3 , so that the 



characteristic mass given by eq. (40) becomes ~ 100 Mq. The evolution at higher 
densities is similar to what described before. The main difference is the predicted 
lower masses that would result from the fragmentation of primordial gas into 
smaller clouds of typical solar /subsolar mass (e.g., Nakamura & Umemura 2002; 
Nagakura &: Omukai 2005; Johnson &: Bromm 2006). 

6 INTERACTIONS OF MOLECULES WITH THE CMB 

The study of the interaction of primordial molecules with the CMB is important 
in view of the possibility to observationally detect spectral distortions or spatial 
anisotropics imprinted by molecular transitions on CMB photons as they travel 
from the last scattering surface to a present-day observer. Before the COBE 
mission (Smoot et al. 1992), much of the interest about primordial molecules as a 
source of opacity was motivated by the lack of evidence for spatial anisotropics in 
the CMB, leaving the possibility that an absorbing layer of molecules extending 
over a significant range of redshifts could smear out any primary anisotropy 
produced at higher redshift and possibly generate detectable secondary spa- 
tial/spectral features (Dubrovich 1977,1994; Maoli, Melchiorri & Tosti 1994; 
Maoli et al. 1996; Puy & Signore 2002). 

Starting from the late 70s, the main effect investigated was resonant scatter- 
ing by molecules in protogalactic clouds moving at some peculiar velocity with 
respect to the Hubble flow (Dubrovich 1977), a process analogous to the Thom- 
son scattering of CMB photons by free electrons (Sunyaev & Zel'dovich 1970). 
Resonant scattering by molecules, like the Sunyaev-Zel'dovich effect, is an elas- 
tic process in which a low energy photon [hv <C m e c 2 ) is absorbed and then 
re-emitted at the same frequency but in a different direction, isotropically in the 
rest frame of the absorber but anisotropically in the observer's frame. If the scat- 
tering particles have a density n s and a peculiar velocity v pec along the line of 
sight (positive away from the observer), the induced spectral/spatial anisotropy 
in the CMB temperature to the lowest order in v pec /c is 

^ = -^£(l-e--) (41) 
T r c 

where 

n s a u d£, (42) 

is the optical depth of the protogalactic cloud. For free electrons, a u is indepen- 
dent of frequency and equal to the Thomson cross section <tt = 
whereas for electrons bound in molecules the scattering is frequency-dependent 
and resonant at the frequency of a molecular transition: for a single transition 
with frequency v u [ between two levels with density populations m, n u ;, statistical 



The Dawn of Chemistry 



35 



weights gi, g u and Einstein coefficient A u \, the absorption cross section is 

; A ul 9 ^(u - u ul ) (43) 



87ri/ 2 V g u nij g\ 

where (j>(y — v n \) is the line profile. The important fact here is that the absorption 
cross section of electrons bound in molecules is several orders of magnitude larger 
than the Thomson cross section of free electrons. Resonant scattering can also 
be greatly enhanced by luminescence, i.e. the absorption of a high energy photon 
corresponding to a rovibrational transition followed by the re-emission of several 
lower energy photons (Dubrovich & Lipovka 1995). Since at z ~ 300-500 the 
frequencies of rotational transitions fall in the Rayleigh-Jeans region of the CMB 
whereas the frequencies of rovibrational transitions fall around the peak, where 
the photon density is maximum, luminescence can result in an increase of AT r /T r 
at low frequencies by orders of magnitudes for molecules possessing low-energy 
rotational transitions like LiH and H2D + (Dubrovich 1997). 

The limits of integrations in equation ( |42[ ) depend on the size of the proto- 
galactic cloud where molecules are located. In the case of a uniform distribution 
of molecules in the primordial gas, with the abundance depending on redshift as 
described in the previous sections, is convenient to transform the integration over 
space in an integration over redshift, dJL = cdz/H(z), where H(z) is the Hubble 



parameter defined by equation (21). In practice, the absorption (or emission) 
of photons observed today with frequency v occurs in a small range of redshifts 
Azi « (1 + Zi)Avth/v u \ around the interaction redshift z\ defined by the condi- 
tion z/(l + Zi) = z^ui. Since spectral lines expected to arise in the dark ages are 
narrow (Afth/f u i *C 1), the region contributing to the optical depth is also small 
(Azi/z; <C 1) and can be approximated by a Dirac delta in evaluating the integral 
in equation p2|) (for details, see e.g. Bougleux k, Galli 1997). 



6.1 Optical Depth 

In a cold gas made mainly of H and He very few atomic transitions can be excited 
at low temperatures. A notable exception is the hyperfine 21 cm transition of 
neutral hydrogen, with an energy difference between the F = 1 and F = levels 
corresponding to 0.068 K. This spin-flip transition is expected to have resonantly 
absorbed CMB photons in the post-recombination era, producing a broad ab- 
sorption feature at frequencies of ~ 10 MHz (Varshalovich & Khersonskii 1977, 
see also Field 1959). However, the detection of this feature is severely hindered, 
at least from the ground, by ionospheric distortion and reflection of radio waves 
and by the Galactic synchrotron emission that dominates the CMB by several 
orders of magnitude in this frequency range (Jester & Falcke 2009). The sensi- 
tivity of this transition to a variety of processes that can affect the temperature 
evolution of the gas before and during reionization, when the 21 cm signal (either 
in absorption or in emission) is redshifted to more convenient frequency bands 
around ~ 100 MHz, has made this line a prime target for the new generation 
of radio interferometers like LOFAR and SKA (for a recent review of 21 cm line 
cosmology, see Pritchard & Loeb 2012). 

In addition to the 21 cm line of neutral H, the 6708 A resonance transition 
between the ground state 2 S and the first excited state 2 P of neutral Li can also 
provide a source of opacity of the primordial gas, as first proposed by Loeb (2001). 
Stancil et al. (2002) found that this process could result in a possibly detectable 
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suppression of the power spectrum amplitude of the CMB. However they did 
not consider the dramatic impact of non-thermal photons on the recombination 
history of Li, as discussed in Section lithiumrecombination. When this effect is 
taken into account, the abundance of neutral Li decreases by orders of magnitude, 
and the optical depth of the Li 6708 A line is reduced to only a few 10 -5 (Switzer 
& Hirata 2005). 

Unlike atoms, molecules possess rovibrational transitions that can be excited 
at the low temperatures characteristic of the pregalactic gas. Early calculations 
of the contribution of primordial molecules to the opacity of the universe were 
heavily affected by large uncertainties on the fractional abundance of the species 
considered. Dubrovich (1977) computed the optical depth of the lowest rota- 
tional transitions of HD + , LiH, para- and ortho-H2D + and found that the latter 
could reach unity at cm wavelengths if H2D + /H ~ 10~ 6 . Special attention was 
given to LiH, because of its high dipole moment (5.89 debyes, Wharton, Gold & 
Klemperer 1960), and therefore its strong rotational and rovibrational transitions 
(see, e.g., Zemke & Stwalley 1980; Bellini et al. 1994; Gianturco et al. 1996). In 
addition, because of the lack of accurate quantal rates for the radiative associa- 
tion of LiH, during the 1980s and early 1990s the conversion of Li into LiH was 



expected to be almost complete, as described in Section 4.5 



Recently, Black (2006) has computed the opacity associated with bound-free 
photodetachment transitions of the negative hydrogen ion H~ on the CMB, and 
found an optical depth of more than 10~ 5 at mm wavelengths, which would 
have observable effects on the CMB. Following Black's suggestion, Schleicher 
et al. (2008) analyzed both the bound-free and the free- free absorption of H , 
including the corrections for stimulated and spontaneous emission, and found 
that the optical depth is dominated by the free-free process for frequencies below 
~ 100 GHz, with a v~ 2 frequency dependence, and by the bound-free process 
at higher frequencies. The value of the optical depth however is only ~ 1CP 9 at 
mm wavelengths, much lower than the estimate by Black (2006). In any case, 
the negative hydrogen ion remains the main source of opacity in the dark ages at 
frequencies below ~ 30 GHz, whereas at higher frequencies the optical depth is 
dominated by the contribution of rotational transitions of HeH + and the 2 S- 2 P 
transition of Li. 

Figure 7 shows the optical depth as function of frequency of H (21 cm line), Li 
(6708 A line), H~, HeH + , HD + , LiH and H2D + , computed with the abundances 
shown in Figures 3a-d. From lower to higher frequencies, the optical depth is 
dominated by pure rotational transitions and rovibrational transitions (free-free 
and bound-free transitions, respectively, in the case of H~). For LiH and H2D + 
only the contribution of pure rotational transitions is present in this range of 
frequencies. 



6.2 Spectral Distortions 

Besides the smallness of the optical depth of the main molecular species, another 
limitation to the detection of spectral distortions in the CMB induced by resonant 
molecular scattering is the low value of the peculiar velocity v pec with respect to 
the sound speed. The magnitude of the peculiar velocity has a maximum value 
of ~ 900 km s _1 for protogalactic clouds of masses 10 9 -10 11 M and is lower 
for larger or smaller masses (see e.g. Maoli et al. 1996). It also decreases with 
redshift as (1 + z) -1 / 2 . Thus, v pec /c < 10~ 3 /(1 + z) 1 / 2 . Contributions to v pec can 
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Figure 7: Redshift-integrated optical depth for several primordial atoms, 
molecules and ions as function of frequency/wavelength. For clarity, the opti- 
cal depth due to the 21 cm line of neutral hydrogen is shown on a separate panel, 
and the contribution of HD + has been displaced upwards by a factor of 10. For 
comparison, the intensity of the CMB and the Cosmic Infrared Backround (CIB) 
are also shown in arbitrary units. Details of the calculations for each species 
can be found in Palla, Galli & Silk (1995), Bougleux & Galli (1997), Palla & 
Galli (2000), Switzer & Hirata (2005), Schleicher et al. (2008). 



also come from internal hydrodynamical flows, like gravitational collapse. The 
maximum collapse velocity is in this case of the order of the free-fall speed, or 
about 100-300 km s~ 1 in the same mass range. Finally, the observation of spectral 
distortions induced by resonant scattering is limited to regions of relatively low 
densities: above some critical density line emission will start to dominate scat- 
tering, erasing any decrement of the CMB brightness temperature produced by 
scattering protoclouds (Basu 2007). For example, densities of about 0.1-1 cm -3 
are sufficient to erase the scattering signature produced by the lowest rotational 
transitions of HD and HD + . 

The intensity and angular size of emission/absorption lines produced by various 
molecules during the evolution of primordial density perturbations with masses 
ranging from 10 9 to 10 12 M were computed by Maoli, Melchiorri & Tosti (1994), 
Maoli et al. (1996). They found that the very initial phases of nonlinear collapse 
provide the best conditions for the detection of resonant scattering. In fact, at 
the moment when the cloud collapse is exactly compensated by the cosmological 
expansion (turnaround), the line width reaches the minimum (i.e. thermal) value, 
and the absorption cross section is the highest (Zel'dovich 1978). The signal 
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thus changes from a relatively low-intensity broad line in the linear phase to 
an intermediate-intensity line of growing broadening during nonlinear collapse, 
passing through a high-intensity thin peak at turnaround. 

These studies prompted a few attempts to detect spectral signatures of pri- 
mordial molecules. The first one, performed at mm wavelengths with the IRAM 
radiotelescope, consisted in a search for redshifted rovibrational transitions of 
LiH emitted by protogalaxies at z ~ 100-200 (de Bernardis et al. 1993). These 
observations gave upper limits of ~ 20 mK on the intensity of LiH lines, only 
slightly above the expected intensity of LiH for optical depths of order unity, 
a value that was later found largely overestimated. Using a chemical network 
based on the reaction rates for Li chemistry known at the time, Palla, Galli & 
Silk (1995) found LiH abundances of the order of 10~ 14 -10~ 13 and a value of 
only ~ 10~ 5 for the optical depth of LiH at mm wavelengths. Bellini et al. (1994) 
suggested to extend the search at cm wavelengths, where the optical depth of 
redshifted pure rotational transitions of the v = level would be about two or- 
ders of magnitude larger. However, no results were obtained by Gosachinskij et 
al. (2002) in their search for spectral signatures of LiH at a wavelength of 6.2 cm 
using the RATAN-600 radio telescope. Using more reliable abundances, the op- 
tical depth of LiH was found to be of the order of only ~ 10~ n and ~ 10 -13 at 
cm and mm wavelengths, respectively (Bougleux & Galli 1997). The search for 
HD lines appears more promising, as the optical depths of the three lowest pure 
rotational lines are much higher than for LiH, although their intensities are still 
below the sensitivities of current facilities like GMT/LMT, CARMA and ALMA 
(Nunez-Lopez, Lipovka & Avila- Reese 2006). Recently, Persson et al. (2010) have 
reported the results of broad band, blind line surveys with the Odin satellite to 
search for the lowest rotational transitions of H2, as well as higher transitions of 
HD and HeH + redshifted from z ~ 20-30, the epoch of the first structure forma- 
tion. Although no lines were detected, the search set some interesting limits on 
the amplitudes of the possible resonant lines and can serve as a pilot study for 
current and future radiotelescopes. 



6.3 Spatial Anisotropies 

In addition to spectral features, resonant scattering from primordial molecules 
can also produce secondary spatial anisotropies in the CMB, leading to changes 
in the expansion coefficients C\ of the CMB power spectrum 

^ « -2r„ (44) 

for t u <C 1 and I » 1 (Maoli, Melchiorri & Tosti 1994; Basu, Hernandez- 
Monteagudo & Sunyaev 2004; Basu 2007). Thus, a promising way to detect 
secondary anisotropies due to resonant scattering is to compare the expansion 
coefficients of the angular power spectrum C\ at different wavelengths to detect 
signals of ~ 0.1 fxK or less, exploiting the sensitivity of satellites like Planck, and 
other ground-based or balloon-borne instruments like the Atacama Cosmological 
Telescope (ACT), Atacama Pathfinder Experiment Telescope (APEX) and the 
South Pole Telescope (SPT). However, since the angular scales of the expected 
anisotropies are small, of the order of a few arcmin at best (Maoli et al. 1996), 
only multipole coefficients with I > 1000 will be affected, making the detection 
difficult. A possible experimental setup for the Multi-Beam Solar Radio Telescope 
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(MSRT, Tuorla, Finland) is described in Dubrovich, Bajkova & Khaikin (2008). 
6.4 Emission by Molecular Formation 

In analogy with H and He recombination, molecule formation also occurs in ex- 
cited states that then decay to the fundamental level emitting photons that pro- 
duce an excess over the pure Planck distribution of the CMB. For example, H2 
is formed by associative detachment in vibrationally excited states that quickly 



decay to the ground state (see, Section 4.2). Since the maximum production of 



H2 occurs at redshifts z < 100, the emission of rovibrational transitions with 
A « 2 fim is redshifted today at A < 200 /xm in the Wien part of the CMB. Co- 
incidentally, this is the same wavelength range where photons produced by Ly-a 
and two-photon decay of H and He during recombination are also redshifted (Sec- 
tion 3.3). The H2 emission has been computed by Coppola et al. (2012, see also 
Khersonskii 1982 and Shchekinov & Entel 1984) for multiquantum transitions up 
to Av = 4. As in the case of H and He recombination photons, a direct detection 
of the H2 formation emission in this wavelength range appears challenging, ow- 
ing to the presence of an infrared background (both Galactic and extragalactic) 
several orders of magnitude brighter than the CMB (see e.g. discussion in Wong, 
Seager & Scott 2006). However, the excess photons produced by H2 formation 
could significantly affect a number of photodestruction rates at high redshift. 

6.5 Absorption by Thermal Nonequilibrium 

At redshift below z ~ 500 matter and radiation are no longer efficiently coupled 
by Compton scattering, and the temperatures of baryon and photons decrease 
approximately according to the adiabatic laws T r oc (1 + z) -1 and T g oc (l + z)~ 2 , 



respectively (see Section 3.3). In these circumstances, molecules act as heat 
pumps, transferring energy from the hotter radiation to the colder gas (Kherson- 
skii 1986; Puy et al. 1993). The effect on the CMB temperature is of the order 



where T ex is the excitation temperature of the transition, defined by the condition 
g\n u /g u n\ = exp(—hv u i/kT ex ). In a collapsing cloud, when the density becomes 
larger than a critical density for thermal equilibrium with collisions, T ex becomes 
equal to T g , and the latter rapidly increases above T r . In this case line emission 
start to dominate (e.g., Mizusawa et al. 2005, Basu 2007). On the other hand, 
in the expanding primordial gas after recombination, when radiative effects still 
largely dominate over collisional excitation/de-excitation, the excitation temper- 
ature of a molecular transition is very close to the temperature of the CMB, 
T g < T ex < T T . This explains the limited absorbing power of a layer of primordial 
molecules and the major obstacle for their detection: those species with high 
dipole moments, and thus high absorption cross sections, have T ex ~ T r , while 
species with low dipole moments have t u <C 1. Puy et al. (1993) considered 
only pure rotational transitions in the v = level of H2, HD and LiH and found 
negligible effects of molecular absorption on the CMB: AT T /T T < 10~ 9 at v ~ 40- 
80 GHz for rotational transitions of H2 and even less for rotational transitions of 
HD and LiH. Also in the case of H~, the relative change in the CMB temperature 
due to the bound-free and free- free processes discussed above is small, because the 
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absorption and the emission processes are close to equilibrium, making the net 
effect of the order of only ~ 10~ 10 over a frequency range from 100 to 1000 GHz 
(Schleicher et al. 2008). As already mentioned in Section 6.1, the 21 cm hyperfine 
transition of neutral hydrogen represents an exception to this behavior. Because 
of its extremely small Einstein coefficient (^4io = 2.85 x 10 -15 s^ 1 ), collisions are 
able to couple T ex (called in this case "spin temperature" ) to T g down to z « 100, 
when the Hubble time becomes longer than the lifetime of the F = 1 level and 
T ex approaches T r until reionization occurs (Purcell & Field 1956; Pritchard &: 
Loeb 2008). 



7 SUMMARY POINTS 



• Assuming that cosmology and standard Big Bang nucleosynthesis have en- 
tered the precision era, the dawn of chemistry can now be accurately fol- 
lowed from the recombination epoch to the end of the dark ages. Chemical 
networks containing more than ~ 250 reactions for ~ 30 species have been 
developed to make quantitative predictions on their abundances. Most of 
the remaining uncertainty can be attributed to the partial knowledge of 
cross sections and rate coefficients at the relevant energies and tempera- 
tures. 

• Due to the lack of dust grains, primordial chemistry is driven by slow gas- 
phase reactions promoted by collisions with electrons and protons left over 
from recombination at a level of one part in ~ 10 4 . As a result, only simple 
diatomic and triatomic molecules and molecular ions containing H, D, He, 
and Li are formed in trace amounts. 

• The low density of the expanding medium and the harsh radiation field due 
to the cosmic background limit the maximum production of each species. 
While H2 and HD reach freeze-out abundances of < 10 -6 and > 10~ 10 , 
respectively, at z ~ 40, all other molecules/ions have abundances much 
smaller than 10~ 10 . 

• The evolution of H2, H^ , Hjj" and their deuterated isotopologues depends to 
a large extent on the additional contribution to the cosmic radiation field 
of photons produced by the recombination of H and He. These photons 
also reduce the abundance of H~ and Li: the former limits the freeze-out 
fraction of H2, the latter never recombines. 

• Primordial molecules allow the gas to cool, contract, and fragment. While 
the answer to fundamental questions regarding the formation of the first 
structures and stars (fragmentation scale, mass spectrum, maximum mass, 
etc.) does not depend on chemistry alone, a full understanding of the ther- 
modynamical properties of the baryonic matter is a fundamental ingredient 
in all numerical simulations addressing these issues. It is rewarding that 
detailed chemical networks and accurate cooling functions are now integral 
parts of state-of-the-art simulations. 

• The interaction of the CMB with the primordial gas is probably the best 
probe to test the predictions of primordial chemistry and provides a unique 
way to peer into the dark ages prior to structure formation. The effects of 
absorption, emission and resonant scattering on the temperature and power 
spectrum of the CMB have been fully explored theoretically. 
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• The most promising results are obtained for the negative hydrogen ion 
H and the HeH + molecular ion. The free-free process of H~ leads to a 
frequency-dependent change in the power spectrum of the CMB of the order 
of 10 -7 at 30 GHz. HeH + can efficiently scatter CMB photons and smear 
out primordial fluctuations in the frequency range 30-300 GHz, leading to 
a change in the power spectrum of the order of 10~ 8 . Finally, resonant 
scatter by neutral Li dominates the optical depth of the universe in the 
far-infrared. 

8 FUTURE ISSUES 

• Major progress has been achieved by the inclusion of state-resolved chemical 
reaction rates for H2 and H^ , that are formed out of equilibrium in highly- 
excited vibrational levels and then decay radiatively to the fundamental 
state. It is desirable to extend this approach to other species like HD, HD + 
and HeH + . Further improvement comes from recent fully-quantal ab initio 
calculations of reaction cross sections for He and Li chemistry specifically 
performed for applications to the early universe. Extending these calcu- 
lations to other species/processes is an obvious need. It is necessary to 
complement these theoretical studies with laboratory experiments, espe- 
cially at sub-eV energies. 

• While excellent agreement between theory and experiments has been reached 
for some key reactions, most notably the dissociative attachment of H 
leading to H2 formation and the dissociative recombination of H^, a num- 
ber of other important reactions still suffer from either poor knowledge or 
inconsistencies. For example, in the high density phase of collapse within 
dark matter minihalos, the conversion to fully molecular gas by three-body 
reactions is still affected by significant uncertainties. We have provided an 
assessment for this and several other critical reactions. 

• Searches for signals from primordial molecules and for their imprint on 
the CMB are very challenging, in part owing to the extreme weakness of 
the expected lines and also because of the poorly constrained redshift of 
emission. However, the exploration of the dark ages using molecular and 
atomic transitions is already under way with the largest available facilities 
like ALMA and LOFAR, and dedicated satellites like Planck. We can hope 
that with the advent of next-generation telescopes (e.g., JWST, SKA, ELTs, 
etc.) some fundamental breakthrough can be achieved that will shed light 
on the dawn of chemistry and its impact on the birth of the first structures. 
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